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APPENDIX  A 


NTS-2  MASS  PROPERTIES 


The  mass  properties  for  NTS-2,  Tables  A-1  and  A-2,  consist  of  launch  and  orbit 
conditions.  DRY  signifies  NO  hydrazine  in  that  condition.  Hydrazine  usage  can 
be  obtained  from  Table  A-3  and  be  added  to  any  dry  condition  by 

2 2 

I Total  = I+  WTxd  +1^+  x d, 
c c c h h h 

where 

C = Condition 
H = Hydrazine 

d = distance  CG  is  transferred 

Individual  component  mass  properties  are  shown  in  Tables  A-1  and  A-2,  which 
along  with  hydrazine  usage  can  be  added  to  the  last  two  conditions  shown  in 
Figure  A-1. 

EXAMPLE: 

USING  COND  (DRY),  AKM  BURN,  AKM  SEP,  WITHOUT  SOLAR  ARRAY,  GRAVITY  GRADIENT, 
NUTATION  DAMPER,  MAGNETOMETER,  AND  REACTION  WHEELS,  FIND  MASS  PROPERTIES 
WITH  SOLAR  ARRAYS  STOWED  AND  50  PERCENT  HYDRAZINE  USAGE. 


New  C .G. Z . 


EW  X CGZ 
TWT 

(799.7)  (25.69  + (60.6)  (22.00)  + (26) C 31.68) 

886.3 


I 

XX 


I 

yy 


I 

zz 


= 25.61  in. 

70.02  + (799.7)  (0)^  + 11.82  + (60.6)  (0)^  + 0.1  + (26)  (0)^ 


81.94  Slug-Ft. 

_ „ , (799.7)(.08)^  4.  R A.:  4.  (60.6)(3.61)^  4.  r i 4-  (26)(6.07)^ 
DP.jz  + (32.2)(144)  ^ ^ (32.2)(144)  (32. 3)  (144) 

68.45  Slug-Ft.^ 

102  40  + (799.7).(..08)i  + 14  47  + (.60 .6)  (3.61).^  .,3  m 
102.40  + (32.2)(144)  (32.2)(144)  ^ ^ (32.2)  (144) 

120.35  Slug-Ft.^ 
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CENTER  OF  GRAVITY  MOMENT  OF  INERTIA  PRODUCT  OF  I^RTIA 
(In) (slug-ft^) (slug-ft^) 
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Figure  A-1.  Hydrazine  Usage 


APPENDIX  B 


ONE  BODY  MODEL  EQUATIONS 


This  appendix  presents  the  derivation  of  equations  of  motion  used  in  the 
one  body  simulation.  For  Definition  of  Symbols,  see  Figure  B-1. 


Vector  torque  equations  are  given  by 

f rxdF  = / rx  r dm 
dt 


dt 


/ rx  C-^)  dm 


T 


dt 


(B-1) 


where  H is  defined  as 


H = / rx  (-T^)  dm 
at 


X i i 


-Iw  +IIa). 
o o i i ^ 


(B-2) 


where , in  body  0 coordinates  , 

AT  T 

M(a,b)  = a bl  - a b 


-I 


xy 


-I 


-I 


xy 


xz 


xz 


xy 


-I 


xy 


J,  = T ,J,  T,  (T  , is  transformation  from  wheel  1 to  body  o) 
i ol  i lo  ol 
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■7" 


I = f (y^  + z^)din  = I + I , etc. 
X ^ yy  zz  ’ 


I = / xy  dm,  etc. 
xy 


The  rotational  equations  of  motion  are  then  given  by 


^ dtr 


dt 


r I 0)  4"  ^ I . cu . *] 

•-  O o , 1 I-* 


“ '■^o'"o  ^ ^i'^o^  “o^  ^^o“o  ^ ^i'^i^ 

1 1 


T = 1 


(B-3) 


Wheel  Torques  may  be  described  by 


/ x^xdF  = / x^x  — 2 
dt 

d ‘^^i 

“ d?  ^ ^i^  ^ 


dH, 


T.  = 


i dt 


Wheel  momentum  is  given  by 


^ ^ ^i^  -ZF  ^i 


= / x^x  X Xj^] 


dm 


< 

I 


Xix[i»)i  X x^  + X x^]  dm 


J U)  + J 0) 
11  1 o 


(B-4) 


Differentiating 


T.  » -1  <0.  + J.i)  + u)  X J.ii).  + 0)  X J.w 
1 ll  lo  o 11  o lo 


(B-5) 


Torques  about  the  wheel  axis  are  given  by 


^1  • •'l^  ^ •'l-o 


(B-6) 


- J^a>^ 


(B-7) 


where  Is  the  sum  of  motor,  bearing  and  windage  torques  on  wheel  1. 
Equations  B-7  and  B-3  are  combined  to  compute  the  body  rate  derivatives. 


Id)  + rT.  + (jo  xH 
0 0^1  o 


(B-8) 


Kinematics 


The  Inertial  coordinate  system  and  the  Euler  sequence  are  consistent  with  the 
Rockwell  convention  presented  in  the  AVCS  splinter  meeting  notes  at  the 
preliminary  design  review.  Specifically,  Rockwell  specified  a yaw  pitch 
roll  Euler  sequence  with  the  equations: 

6 = w cos  + q cos  i})  - r sin  cfi 

0 = (o)  sin  Jp  sin  6 + q sin  4i  + r cos  <}>)  (1/cos  0) 

= p + iji  sin  0 + 0)  sin  <p  cos  0 
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whe  re 

u = orbit  rate 
P»q»r  = body  rates 
<(),0,ip  = Euler  angles 
4),6,ii)  = Euler  rates 

These  equations  are  consistent  with: 

1.  A local  vertical  system  with  z axis  pointing  toward  the  earth  and  x 
axis  in  the  orbit  plane  and  positive  on  the  direction  of  motion. 

2.  Euler  angles  specifying  body  orientation  relative  to  the  local  vertical 
frame . 

The  inertial  coordinate  system  is  the  Geocentric  Equatorial  Frame  (e.g.  , i 
axis  positive  North  and  x axis  positive  in  the  direction  of  the  Vernal 
Equinox).  The  local  vertical  system  is  consistent  with  Rockwell's  (e.g., 
z axis  positive  toward  the  earth  and  x axis  in  the  orbit  plane). 

Two  sets  of  Euler  angles,  both  having  a yaw  pitch  roll  sequence,  are  used. 

The  first  set  PSI , THETAl,  PHIl  define  the  inertial  to  body  transformation  E. 
The  second  set,  PSILV,  THETA,  PHI  specify  the  orientation  of  the  body 
relative  to  the  local  vertical  frame  and  should  be  consistent  with  Rockwell's 
Euler  angles.  The  inertial  Euler  angles  are  integrated  in  the  simulation 
via  quaternions. 

Linear  Motion 

The  simulation  uses  two  subroutines f POSITN  and  RADVEC)  to  compute  the 
I ephemerides  for  the  satellite  and  the  sun.  The  ephemerides  for  the  satellite 

and  the  sun  are  calculated  using  the  six  Kepler  elements:  semi  major  axis  (a), 

i 

I 

I 


eccentricity  (e) , inclination  (i) , right  ascension  of  the  ascending  node  (n) , 
argument  of  perigee  (u) , and  mean  anomaly  (M) . The  update  is  a simple  rate 
term  on  the  latter  three.  Table  B-1  presents  the  six  elements,  their  units 
and  values  for  the  vernal  equinox  condition.  The  rate  terms  are  derived 
in  the  program  from  the  fixed  parameters  and  are  given  by 

fl  ^ v/a  - ^ (1.5  sin^  i-1) 

[a(l-e  ) ] 


n = 


kV  u/a^ 

laa-ehf 


cos  1 


[a(l-e2)]2 


(2.5  sin^  i-2) 


where 

p = Gm  = 626  30.3949 

2 

k = (first  term  of  perturbing  force  = -19255.96124 

due  to  Earth's  oblateness) 


The  rate  terms  for  the  sun  are  fixed  and  are  included  in  Table  B-1.  The 
simulation  updates  M,  u by 

M(t)  = M(t  ) + M • (l-t  ) - TT  < M < TT 

o o — 

0(t)  = 0(t  ) + • (t-t  ) - TT  < 0 < TI 

o o — 

<;l)(t)  = (u(t)+to’  (t-t)  -Tr<u)<TT 

o o — 

A seventh  term,  the  Greenwich  Hour  Angle  (GHA) , which  relates  the  earth- 
fixed  coordinate  system  to  the  inertial  system,  given  the  six  elements  is 
updated  via: 

GHA(t)  = GHA(t  ) + u ■ (t-t  ) 
o e o 

0)  = earth  rotation  rate 

e 
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TABLE  B-l.  GPS  AND  SOLAR  EPHEMERIDES 


STANDARD  VALUE 

SYMBOL 

NAME 

UNITS 

GPS 

SUN 

a 

Semi  Major  Axis 

Miles 

14342. 

1* 

e 

Eccentricity 

— 

0. 

.0167 

i 

Inclination 

Degrees 

63. 

23.443 

M 

Mean  Anomaly 

Degrees 

0. 

88.068 

U 

Right  ascension 

of  ascending 

node 

Degrees 

0. 

0. 

(l) 

Argument  of 

pe  ri  gee 

Degrees 

0. 

270. 

M 

Degrees /sec 

— 

1.14074x10"^ 

h 

1 

Degrees /sec 

— 

.544773x10“^ 

0) 

Degrees  /sec 

— 

0.0 

OIA 

Greenwich  Hour 

Angle 

j 

Degrees 

0.0 

A 

Astronomical  unit 


APPENDIX  C 


THREE  BODY  MODEL 


1 

I 

I 

I 

I 

1 

[ 

I 

I 


STRUCTURE 

The  overall  "three"  body  structure  is  as  shewn  in  Figure  C-1.  It  differs 

from  that  given  in  Reference  C-1  in  the  addition  of  the  gravity  gradient  rods 

and  tip  mass  components.  The  vector  r^^  extends  from  the  center  of  the 

space  vehicle,  to  the  support  point  of  rod  k.  The  vector  Pj^  extends 

from  the  support  point  to  the  center  of  tip  mass  k.  The  vector  Xj^  extends 

from  the  center  of  tip  mass  k to  a differential  mass  element.  Since  all 

torques  will  be  referenced  to  the  total  system  center  of  mass,  CM,  the 

vector  r will  contain  the  effects  of  center  of  mass  variation  due  to  rod/tip 
o 

mass  movement. 


EQUATION  DERIVATION 

In  modeling  the  dynamics  associated  with  the  gravity  gradient  rods  and  tip 
masses,  the  form,  procedure,  and  notation  developed  in  Reference  C-1  was 
followed  as  closely  as  possible. 


The  system  momentum,  with  respect  to  Inertial  space  is  given  by: 

H = I 0)  + E.I.o).  + I.I.u),  + I,  H,  + E H 
oo  ill  jjj  kk  nn 

where  the  inertia  tensors 


(C-1) 


lo  = ^o  ^I'^i  ^k\  ^n'^n  ~ ^ Ej^m^M(r^  ,r^) 

+ ^j“jM(rj+pj,r.+Pj)  + Z^mj^M(rj^+p^.r^+p^)  + E^m^M(r^+p^  ,r^+p^) 

(C-2) 


II  = 


(C-3) 
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GRAVITY  GRADIENT  ROD  SUPPORT  POINT  P,^ 
SOLAR  PANEL  PIVOT  POINTS  P^ 

BODY  CM's:  CM.,  CMj,  CM,^;  CM  = TOTAL  CM 
MASSES:  , nij,  m|^,  m = TOTAL  MASS 

INERTIA  TENSORS  (ABOUT  CM's)  J.,  J j , 


Figure  C-1.  Three  Body  Structure 


( 


. Jj  + Y(r 


(C-4) 


^li 

M(a,b)  - b^al  - ba^ 
and  where  is  conputed  from 

+ Z iii(r+r  + p ) = 0 (C-5) 

n n o n n 

If  we  assume  a massless  rod  and  the  sphere  is  a point  mass,  then 

J = 0 (C-6) 

k 

The  momentum  may  be  represented  by 

«Pk  - - - 

«k  = V’^k  -^k^  ^ ^ ^k^’^o'^\^k‘'V  ^ rr 

Since  we  have  assumed  the  tip  mass  may  be  represented  by  a point  mass,  the 
Second  term  in  equation  c-7  tnay  be  set  to  zero. 


The  torque  equations  which  are  the  ultimate  goal  of  the  model  are  derived 
from 


T 


dH  6H  - 
dE  " ^ “o 


X H 


(C-8) 


Or,  differentiating  (1) 


I.iL  + f 


o o 


u + I.I.U,  + T.I.w,  + 1,  H,  + EH  + W xii  (c_9) 
o iii  jjj  kk  nn  o ^ / 


where  the  Inertia  tensor  derivatives 

+ £jmj(M(rj+Pj  , ;Pj)) 

k !:„»„(M(r„kp„,.'„)  k M^r„kp_^.i„))  (C-10) 
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where 


ij  - jj  + j - o;'  + '“j> 


(c-ll) 


Pj  = (.J  X Pj 


Wk  + ^n'"n'^n) 


(C-12) 

(C-13) 


The  derivatives  of  the  rod/tip  mass  momentum  shewn  in  equation  c-9  is  given 
by 


\ = “k 


6t  6 1 


5Pi 


6 Pv 


WVf^k^  *-2 


6t 


(C-14) 


6 Pi 


— ^ is  the  acceleration  of  tip  mass  k with  respect  to  the  space  vehicle. 
6t 

At  this  point  we  introduce  rod  flexure  considerations.  Consider  the  location 
of  the  tip  mass  to  be  defined 


‘'o  \ 


(C-15) 


where 


Pj^  = initial  position  of  tip  mass  k 
o 

= displacement  of  tip  mass  k due  to  flexure 
Differentiating  (15)  with  respect  to  body 


j-r  5n 

6P,  i 


5t  5 1 


and 


5^P, 


5t 


6t" 


(C-16) 


(C-17) 
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( 
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The  flexure  model  for  the  gravity  gradient  rods  is  described  in  Appendix  I. 
It  is  a 2-node  model  which  may  be  described  by 


(C-18) 


Cj^  is  the  damping  matrix  given  by 


where  Cj^  and  C,  have  been  defined  in  Appendix  I as 
a D 

201.06x10"^  0 -.75x10"^  0 

0 201.06x10“^  0 -.75x10 

-.75x10"^  0 12.17x10“^  0 

0 -.75x10“^  0 12.17x10 
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and  Cj^'  is  a scalar  used  for  a conversion  factor  and  for  changing  dancing 
ratios.  NOTE:  It  has  been  assumed  that  there  are  no  deflections  in  the  Z 
direction,  thus  there  are  only  two  degrees  of  freedom  at  each  node. 


K,  is  the  stiffness  matrix 


where 


is  a conversion  factor  and 


46085,^ 

El 

y 

1728^ 

JG 

(■%  -■’x)’ 

0 

1440^^ 

El 

y 

0 

0 

4608^^  . 

-2.  1 

2 

0 

El  ^ 

X 

JG 

o 

El 

X J 

1440)?.^ 

El 


y 


0 


1440«: ' 
El 


0 


5 76g^ 
El 

X -< 


= rod  half  length  = 30  ft  (fully  deployed) 


Ely  = rod  bending  stiffness  about  y axis 


3000110  percent  lb-in 


El  = rod  bending  stiffness  about  x axis  3000+10  percent  lb-in 

X 2 

JF  = rod  torsional  stiffness  25  lb-in 

X > 5^0  “ initial  x and  y components  of  ' representing  unstressed 

position  from  x,  y axes  respectively,  liue  to  thermal  and 
manufacturing  curva  tures . 

~P' 

q is  the  angular  acceleration  of  the  spacecraft  which  is 
.f  J translated  to  the  node  point  of  rod  k by  transformation  Tj^. 
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+ 2^sgn(r^ 


■|jk^  + 2Asgn(rj;^)J 
0 


[■ 


\ + J,sgn(r^  ) 

z z 


(y  - n ) 

^ o y 

\ \ 

^’^k  \ ) 

\ \ 

(y  - n ) 

k k 

\ - X ) 

X D,  b, 


where 


r ,'  = r + r, 
k o k 


The  external  forces,  > ^re  composed  of  gravity  gradient  forces  and 

solar  pressure  forces.  The  solar  pressure  forces  acting  on  the  tip  masses 
are  defined  in  Appendix  G. 

The  gravity  gradient  forces  are  discussed  belcw. 

TIP  MASS  1 
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In  a central  inverse  square  gravitational  field  the  force  acting  on  a particle 
of  mass  m is  given  by 


-mkR 


(C-19) 


where : 

m = mass  of  particle 
k = gravitational  constant 

R = vector  from  the  particle  to  the  nass  center  of  the  attracting  body 


If  we  assume  two  tip  masses,  the  vector  R would  break  dcwn  into  R^  and  R2  as 
shown  in  Figure  C-2. 


R 


1 


R + r: 

o 1 


(C-20) 


R 


2 


4-  R2 


(C-21: 


where  R'  and  R'  are  the  vectors  from  the  center  of  mass  of  the  spacecraf t/rod 
System  to  the  center  of  mass  of  the  tip  masses  and  R^  is  a vector  from  the 
center  of  the  earth  to  the  center  of  mass  oftthe  space  vehicle.  Referring  to 
Figure  C-1 


+ 

+ 


Thus 


-m^kf 


-m^kRj^ 


-m2kR2 


(C-22) 

(C-23) 

(C-24) 

(C-25) 
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Now 


= (R+r+r,  +p,  )’(R+r+r,  +p,  ) 

o o ° ° 

° I ^ I ’'o  ^k  I ^ ^k  ^k  ^ 


= R 


(1  + 2R^  . (r^  + + p^^^)  (?^  + + p-^^^)  ) 


R 


R 


The  last  term  may  be  neglected  leaving 


(1+  2R^  • Cr^-V  r^.  +P,  )) 

I®,!  = IR  I 5 — ' — 

' 1'  ' o'  IR  |2 


Take  the  -3/2  power  of  both  sides 

(l  + 2»<,  • rr„+\  +5 
^7- 


Using  the  binomial  theorem 

(1 


3 R^  • (r^  + r^^  + pj^^)  ) 


IR 


Therefore 


m/k  - 3 Ro  ■ ^’^o  + ’^k,  P'k  )) 

3 ‘^‘^o  ’'o  \ Pk  ^ 2 

1R_I  ° ° S iR  |2 


and 


m k ^ ’^o  ’ ^’^o  ’^k  Pk 

— ^ (S  + I + + J,  ) 3 — 

|R  |3  ° ° '‘2  >‘2  |R  |2 

I o'  'o' 
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The  forces  F.  and  F.  can  be  separated  into  a gravity  component  and  a 
G 

gravity  gradient  component  or 


^ ^ + F^ 

G G gravity  G gravity  gradient 


where 


and 


m,  k R 

? = 1 o 

^G  gravity  ]r  I ^ 

' o ' 


G (gravity  gradient) 


, (l-3R*(r+r,  +p,  ) 

m^k  _ ' o o k ^k/ 

— ^(r  + r + p,  ) T ^ 

|R  1^  ° |R  1^ 


R (3R  ’(r+r,  +p,  ) 

o o o k^ 


ml,  -R  (3R  • (r  + r,  + 0,  )) 

m^k  _ oook  lv 

3 ^^o  \ '^k  ^ 2 — 

|R  I 1 1 |R  I 

'o'  'o' 


Over  the  n body  model  the  gravity  forces  will  cancel  out,  leaving  only  the 
gravity  gradient  forces. 

Summation  of  forces 


,2- 

6 Pj^ 


1 


mi  ^ = F^  + F^ 

. t G s 


(gravity  gradient) 


Where  F^  = solar  force 
s 


All  vectors  are  expressed  in  body  axis,  so  R = T, „ 

^ ■'  ’ o LV  to  B 


0 

0 

R 


\v/B  is  the  transformation  from  local  vertical  axes  to  body. 
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APPENDIX  D 


NUTATION  damper  MODEL 

This  appendix  summarizes  the  nutation  damper  model  implemented  for  the 
digital  simulation  evaluation  of  GPS  stability  and  control  analysis. 

The  nutation  damper  is  used  for  passive  GPS  spin  axis  stabilization  during 
the  initial  mission  phases  by  dissipating  the  spacecraft  precession  energy. 
The  nutation  damper  adopted  for  GPS  application  is  of  the  type  of  ball  in 
gas  filled  tube.  The  tube  is  made  of  aluminum  and  the  ball  of  tungsten 

Carbide.  The  gas  is  a 9 : 1 ratio  mixture  of  nitrogen  and  helium  at  1 atm. 
pressure.  Physical  parameters  of  this  daiipef  is  shown  in  Figure  D-1. 

A similar  damper  of  this  type  was  used  in  the  Telstar  satellite.  Reference 
D-1.  Discussions  on  the  design  and  analysis  of  ball-in-tube  type  nutation 
dampers  can  be  found  in  Reference  D-2. 

Derivation  of  the  torque  equations  presented  here  follows  the  same  procedure 
of  Reference  D-2  which  contains  numerous  typographical  errors.  The  results, 
in  the  torque  form,  are  compatible  for  integration  into  the  three  body  model 
described  in  Appendix  C. 

DAMPER  MODEL  DERIVATION 

To  facilitate  the  model  derivation,  we  shall  define  a vehicle  fixed  B-frame 
whose  origin  is  at  the  vehicle's  center  of  mass  (excluding  damper  ball)  with 
the  (X„  - Y - Z„)  axes  aligned  with  its  principal  axes.  Also  in  the  fol- 
lowing  consideration  of  vehicle  rotational  motions,  it  is  assumed  that  the 
translational  motion  of  vehicle's  center  of  mass  resulted  from  reactional 

I 
( 

I 
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forces  of  the  danper  ball  can  be  neglected.  By  reference  to  the  geometrical 
description  of  damper  location  in  B-frame,  the  position  of  center  of  ball, 

R,  can  be  written  as,  Figure  D-2, 


0 

-a  + Y cos  a 
b + Y sin  a 


(D-1) 


where 

Y = 159  in.  = 13.25  ft. 

(a,b)  = coordinates  of  tube  center  of  curvature 
Equation  (1)  inplles  that  Y„  - Z„  plane  contains  the  danper  tube. 

— — D 


Let  the  vehicle  moment  of  inertia  be  such  that 

r ^XX  “ Vf  ° ^ 

I 'zz  = c 


The  rotational  energy  of  the  vehicle  becomes, 

Ty  - 4 


CD-3) 


The  velocity  of  damper  ball  can  be  written  through  differentiation  of  (1)  as, 
b^  ’ b^  " ^BI^^  b^ 

with  (R^j^)  =■  Tgj^ 

= ^Bl"  ^ "bi" 

= S (r\)  + T T (r\) 
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Zg  = SPIN  AXIS 

cMv  = VEHICLE  CENTER 
OF  MASS 

DAMPER  CONTAINED  IN 
PLANE  (Yg  - Zg) 


Figure  D-2.  Nutation  Damper  in  Vehicle 
Body  Fixed  Frame 


or 


(D-4) 


In  component  forms,  (4)  becomes, 


v„ 

X 

v„ 

- 

Y 

Z 

0 

"2 

-W„ 


-^z 

0 -W^ 

w„  0 


0 

-y  sin  a ci 
y cos  a a 


X 


0 

—a  + Y cos  a 
b + Y sin  a 


(D-5) 


Hence  the  translational  kinetic  energy  of  the  damper  ball  is 
H M (v/  t + v/) 

= ^ M {[-(-a  + y cos  a)  + (b  + Y sin  nt) 

+ [-(b  + , sin  a)  sin  a cl]^ 

+•  f(-a  + Y cos  a)  + y cos  a Y (D-6) 

where  M • mass  of  the  ball 
4 3 

* 3 'Vi,  p * ^^3  gm  = 0.018  slug 

3 

p » density  of  tmgsten  carbide  (WC)  »=  15.7  gm/cm 


Assuming  the  ball  is  rolling  without  slipping  in  the  tube,  the  magnitude  of 
the  angular  velocity  of  the  ball  is  related  to  the  translational  velocity 
of  the  center  of  the  ball  as 


1^2,2  2^2  2 .2-4 

[y  a 01  Y cos  a a j 

^b 


Ya 

^b 


where 


Yb 


radius  of  the  ball  = 0.625  in. 


(D-7) 
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Hence,  the  rotational  energy  of  the  ball  is 

2 4 (|m 


- i M 2 .2 

- y M Y a 


(D-8) 


The  drag  on  the  ball  is  viscous.  Reference  D-3,  since  the  tube  is  gas  filled, 
the  clearance  between  ball  and  tube  is  small  and  the  ball  moves  much  slower 
than  the  speed  of  sound.  Hence,  the  viscous  drag  on  the  ball  can  be 
represented  by  Rayleigh's  dissipation  function  as 


F - (i\)^  (i\> 

where  is  a constant  to  be  described  later. 


The  total  energy  of  the  system  then  becomes 


(D-9) 


T = T + T,^  + T 
V b Y 


= ^ (y^  + b^  + a^  + 2b  y sin  a -2a  y cos  a)} 

X 

^ (b^  + 2b  Y sin  a + Y^  sin^  a)} 


2 2 2 2 
(he  + (a  - 2a  Y cos  a + y cos  a)  } 

2 

(-ab  - a Y sin  a + b y cos  a + y sin  a cos  a)  } 
2 

+Wx  a {M  (y  - a y cos  a + b y sin  a) } 

2 7 2 

My^} 


(D-10) 


With  the  total  vehicle  energy  given  in  (D-IO)  and  the  damper  dissipation  func- 
tion given  in  (D-9),  the  equation  of  motion  of  the  vehicle  and  the  damper 
system  can  be  obtained  as 
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I 


0 


d_  3T  8T 

dt  ~ Z 2Wy 

d_  _ u iT_ 

dt  3Wy  X 


A-  9T_  _ u 3^ 

dt  aw^ 


i_  _ ill 

dt  3d  ~ 3a 


+ w -— 

Y 2W^ 


+ % - 0 


+ -^  =0 

9oi 


(D-11) 


The  rotational  equation  of  motion  of  the  vehicle  can  be  written  in  the  form 
of  Euler's  rigid  body  equation  as 

AW^  + (C  - A)  W^Wy  = 

AWy  - (C  - A)  W^W^,  = Ty 

AW^  = '^z  (D-12) 

where  the  torques  T^,  Ty  , and  are, 

T = -M  [W  (y^  + b^  + a^  + 2b  Y sin  a - 2a  y cos  a) 

+d  (2Wjj  + d)  (by  cos  a + a y sin  «) 

2 

+a  (y  -ay  cos  a + b y sin  a) 

2 2 2 2 2 2 
^^y'^Z  ~ - 2a  Y cos  a - 2b  Y sin  a - y sin  a + y cos  a) 

2 2 2 
+ (Wy  ~ ^z  ^ + a Y sin  a “ b y cos  a.  - y sin  a cos  a)]  (D-13) 

Ty  = -M  [Wy  (b^  + 2b  Y sin  a + y^  sin^  a) 

+ (W  - W Wy)  (ab  + a Y sin  a - b y cos  a - y^  sin  a cos  a) 

I Z X I 
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t 


+ W W (y^  + + 2b  y sin  a - y^  cos^  a) 

2 

+ 2Wy  a (b  y cos  a + y sin  a cos  a) 

2 2 

+ d (b  Y sin  a + y sin  a) ] 

• 2 2 2 

= -M  [W^  (a  - 2a  Y cos  a + Y cos  a) 

. 2 

+ (W^  + (ab  + a Y sin  a - b y cos  a - y sin  a cos 

2 2 2 
+ (-a  + 2a  Y cos  a - y cos  a) 

2 

+ d (a  Y sin  a - y sin  a cos  a) 

- 2Wy  d (b  y sin  a + y^  sin^  a)] 


(D-14) 


a) 


(D-15) 


With  ball's  equation  of  motion  in  tube  given  as, 

^ ya  + (y  “ a cos  a + b sin  a) 

2 2 

- (b  cos  a + a sin  a)  - W^  (b  sin  a +•  y sin  a cos  a) 
2 

“W^  (aslna-Ysina  cos  a) 

2 2 

- (a  cos  a+b  sin  a-  y cos  a + y sin  a) 


+Wd  (as in  a+b  cos  a - a sin  a - b cos  a) 

.A 

+ (^  ) d = 0 


(D-16) 


The  drag  force  on  the  ball  consists  of  two  components:  the  pressure  drag 
and  the  viscous  drag.  For  the  case  where  the  gap  between  tube  and  ball  is 
small,  as  in  the  GPS  nutation  danper,  and  the  speed  of  the  ball  in  the  tube  is 
small  conpared  with  the  speed  of  sound,  it  has  been  determined  in  Ref.  D-3 
that  the  drag  force  is  primarily  due  to  pressure  drag.  The  drag  force  D, 


1 
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can  be  written  as 


i 

I 


D = 


135  ^ ^ ) 

4x64  \ 5/2  / 

' g ' 


where 

[j  = viscosity  of  gas  in  tube 


g = (y^  - Y^j)y^j  = gap  parameter 
= inside  radius  of  tube 
S = speed  of  ball  relative  to  gas 
Cp  = dissipation  coefficient 

= 135  ^ "^b 

4x64  5/2 

g 


For  GPS  nutation  damper, 

p = 178.1  X 10"^  X 2.089  x lO"^  Ib^  sec/ft^ 
g = (0.633  - 0.625)  / 0.625 
= 0.633  in. 


Hence,  value  of  dissipation  coefficient, 
= .0171  lb^/(ft/sec) 


To  handle  the  discontinuity  at  either  end  of  the  damper  tube,  let 
a = -ea  at  a = + a max 

where 

a max  = (8/159)  rad  =2.9  deg. 

£ =0  for  inelastic  collision 

= 1 for  elastic  collision. 
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The  GPS  requirement  on  half  cone  angle  of  nutation  is  to  be  less  than 
2 deg.  The  discontinuity  situation  is  not  expected  to  be  encountered  often 
if  the  2 deg,  requirement  is  met.  The  value  of  £ = 1 would  give  a 
conservative  result  in  that  less  precession  energy  is  dissipated  at  bottoming. 
Hence,  use  of  £ = 1 is  recommended  for  this  reason. 

NUTATION  DAMPING  TIME  CONSTANT  VERIFICATION 


Nominal  Solution  Undamped  Rigid  Body  Motion 


Consider  a rigid  body  with  symmetric  moment  of  inertia  about  the  traverse 
axes,  i.e.  , 1^^  = Iyy  a torque  free  situation.  The  Euler's  equation 
is  reduced  to: 

= -A 

' "y  ■ * "z“x 

= 0 (D-17) 


whe  re 

A = 


A 


(n 


1) 


The  solution  of  the  above  equation  can  be  written  as: 

W = -B  sin  fit 
X 

< = B cos  fit 

_"z  ■ “z^’ 

with  fi  = A W^CO) 


(D-18) 

(D-19) 
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The  angular  momentum  vector  is  written  as: 


- ^XX^X  -B  '''  ^YY  ^Y  -B  ^ZZ  ^Z  ^ 


(D-20) 


Since  there  is  no  external  torque,  the  angular  momentum  is  constant.  The 
nominal  motion  of  the  rigid  body  can  be  visualized  as  a constant  precession 
about  the  angular  momentum  vector  with  a precession  rate  and  a nutation 
angle  6 which  is  defined  as: 


tan  e = — = 


2+1  2w  2 

X ^YY  \ 


^zz^z 


(D-21) 


Therefore : 


B 


ZZ 

X tan  e I W 
\ ^YY  ' 


(D-22) 


For  small  value  of  S,  (6)  becomes 
I. 


B 


X 6 I W 

i^YY  ' ^ 


(D-23) 


A relationship  between  the  rigid  body's  angular  momentum  and  rotational  energy 
useful  for  later  development  of  the  energy  sink  approximation  is: 


[ 


T = ^ (W/  + + J5l„W, 


2 «t 


H 


XX 

2 2 2 
H = H + H 
t s 


ZZ  Z 


^’’XX  ^^ZZ 


(D-24) 


Hence  H. 


From  (5) 


I I ' 
XX  ZZ 

^ZZ  " ^XX, 


2T  - 


ZZ  / 


(D-25) 


2 «t 

sin^  e = 

H 


^xx^zz 

1 1 

1 21 1_ 

I „ - I 

1 

2 ~ 1 

ZZ  xx 

1 H ZZ 

(D-26') 
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Nutation  Damping  Time  Constant 


The  GPS  nutation  damper  consists  of  a tube  aligned  with  respect  to  the 
vehicle  splr  axis.  A ball  inside  the  tube  is  dissipating  the  nutation 
energy  via  vriscous  drag  force  retarding  the  ball's  notion.  An  approximation 
procedure  kncwn  as  the  energy  sink  concept  assumes  that  the  dissipation  of 
nutation  energy  does  not  change  the  system's  angular  momentum  (which  is 
dominantly  spin  angular  momentum) . Hence  equation  (0-26)  can  be  differentiated 
with  terms  H neglected  in  this  approximation,  i.e., 

21  f 

. ZZ 

2 sin  6 cos  0 f?  = j (0-27) 

(n-l)H 


For  small  value  of  0, 


0 0 


21  T 
ZZ 

(n-l)H' 


A further  approximation  involves 

2 2 2 2 

“ • “s  - 'zz  "z 


Hence 


0 = 


(n-l)IzzWz^0 


(D-28) 


(D-29) 


(D-30) 


The  energy  dissipation  T can  be  expressed  as  the  product  of  the  drag  force 
the  ball  experienced  and  the  velocity  of  the  ball,  i.e., 

T = -fv  (0-31) 


: 1 

34  ' 1 

; i 


with  drag  force. 


f = 


and  the  ball  velocity, 
V = ya 


1 


I 

r 

i 


Hence  the  energy  dissipation  rate , 

f = -Cp  (D-32) 

where  a is  the  position  of  ball  in  tube  defined  as  in  Figure  D-2, 
y = is  radius  of  damper  tube, 

Cp  = is  dissipation  coefficient 


From  analysis  of  nutation  damper,  the  motion  of  the  ball  in  tube  was 
obtained  as 


y ya  + W^(y  ~ a cos  a + b sin  a)  - t'Jy^Cb  cos  a + y sin  a cos  a) 
2 2 

-W^  (b  cos  a + a sin  a)  “ (a  sin  a - y sin  a cos  a) 

2 2 

-WyW^Ca  cos  a + b sin  a - y cos  a + > sin  a) 

+ -jr  a - 0 


(D-33) 


With  the  geometry  of  mounting  such  that  b = 0 and  Introducing  the  small 

2 2 

angle  approximations  sin  a = oi , cos  a = 1 and  neglecting  terms  and  Wy 
as  they  are  small  for  small  nutation  angle  0,  equation  (17)  becomes 
C ^ 

(-|^)a  + (~^)oi  + - a)Qi  = W^W^Ca  - v)  - - a)  (D-34) 

Using  the  expression  of  Wy  as  given  in  equations  (2)  and  (7) 


Wy  ” en  cos  fit 


(D-35) 
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Similarly 


W = -6n  W fi  cos 

A id 


Equation  (18)  can  be  reduced  to  the  following  form: 

PSl  . pu-a;Wj2 

a + 


^1 
7M  I 


- 6 


5 

a + - 


7y 


5(y  - a)  (2n  - 


7 Y 


cos  fit 


(n-36) 


The  solution  of  the  standard  second  order  linear  differential  equation  of 
the  form 


a+2ba+a  a=c6  cos  fit 


has  been  obtained  in  standard  textbooks  as: 


a = 


c6 


V(a^  - + 4b 


— TT  COS  (fit  + (f)  = 6a  cos  (fit  + <^) 
2 „2  o 


Compare  (20)  and  (21) , 


2b  = 5 Cjj/(7M) 
a^  = 5(y  - 


-5(y  - a)  (2n  - nhvs^Kly) 


From  (22)  , 


a = - 9 a n sin  (fit  + ip) 
o 


Sii>stltute  (24)  into  (16) 

f - -6^  Y^  a ^ sin^(fit  +(}.) 
D O 


:d-37) 


(D-38) 


(D-39) 


(D-40) 


(D-41) 
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I 


[• 

I 


I 

I 


{‘ 


r 


The  average  value  over  a full  period,  i.e., 

Qt  = , is 

P 

T = -Ij  Cjj 

1 ? 

since  — / sin  + !f')dt  = h 

P 

Substitution  of  (26)  into  (14)  gives 


The 


2 2 2 2 
o 


l_2(n-l)Iz2Wz  J 

time  constant  of  nutation  damping  can  now  be  written  as: 


where 


° V (a^  - n^)^  + 4b^ 
2b  = 5 Cjj/(7M) 

' a^  = 5(y  - a)Wz^/(7Y) 
c = -(2n  - n^)^^ 

= (n  - l)Wz 


(D-42) 


(D-43) 


(D-44) 


r 
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The  predicted  nutation  damping  tine  constant  for  the  two  cases  of  = 
100  RPM  and  = 10  RPM  are  confuted  as  follows  ; 

Case  1 : 

r = 13.25  ft. 

a « 10.25  ft. 

m = 0.018  slug 

Cjj  = 0.0171  lb^/(ft/sec) 

^XX  “ ^YY  * slug-ft 
I22  • 94  slug-ft^ 

'^Z  = 100  RPM  = 10.47  rad/sec. 

n » 1.22 
Q - 2 . 31 

2b  - 0.6789 

a^  = 17.728 

c^  = -16.84 

a - -1.348 
o 

T - 155.  755  sec. 
c 

Case  2 : 

« 10  RPM  * 1.047  rad/sec. 
n - 1.22 
fl  - 0 . 2 31 
2b  - 0 . 6 786 
a^  - 0.17728 


c 


= -0.1684 


a = -0.84275 
o 

T = 398.497  sec. 
c 

( The  nutation  half  cone  angle  history  corresponding  to  Case  1 and  Case  2 

' considered  above  are  plotted  in  Figures  D-3  and  D-4  respectively.  These 

results  are  obtained  using  the  nutation  damping  simulation  program  where  the 
detailed  dynamics  of  the  damper  ball  and  rigid  body  dynamics  are 
I considered. 

I 


( 

f 

i 
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100  RPM 


'xx  = 


Figure  D-3.  Simulation  Data 


10  RPM 


Figure  D-4.  Simulation  Data 
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APPENDIX  E 


ATTITUDE  CONTROL  ELECTRONICS  SYSTEM  MODEL 


The  unmodified  model  of  the  Attitude  Control  Electronics  system  as 
implemented  in  the  computer  simulation  differs  only  from  that  shewn  on  the 
Spacetac  blueprints  in  the  absence  of  voltage  units  and  vjheel  signal 
monitoring  elements.  It  is  given  in  block  diagram  form  in  Figures  E-1, 
E-2,  E-3,  and  E-4. 


I 

I 

I 

I 

I 

f 

I 
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Earth  Sensor 


Figure  E-1.  ACE  Systems  Block  Diagram 
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I Figure  E-2.  Eheel  Command  Transformation  Logic 
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a = 1.57  rad/sec 


b = 2.09  rad/sec 
slope  = 1.0 

(equivalent  to  .14°  db  In  pitch  and  roll, 
2.58°  db  in  yaw,  for  single  axis  inputs) 


Figurr;  E-3.  Nonlinear  Element 
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APPENDIX  F 


TILTED,  CENTERED  DIPOLE  MAGNETIC  TORQUE  MODEL 


The  instantaneous  magnetic  disturbance  torque  is  T = M x B where 

m — — 

M is  the  spacecraft  effective  dipole  moment 
^ is  the  local  magnetic  induction  or  flux  density. 

The  tilted,  centered  dipole  model  is  used  to  represent  the  flux  density. 

The  flux  density  is  a function  of  latitude  and  longitude.  These  parameters 
were  not  being  calculated  in  the  previous  model.  Hence  a new  coordinate 
system,  the  E frame,  was  defined.  The  E frame  has  its  center  at  the  center 
of  mass  of  the  vehicle  with  the  reference  plane  perpendicular  to  the  earth's 
radius  vector. 


The  X axis  points  east,  the  y points  south,  and  the  Z axis  points  towards 
the  center  of  the  earth. 


The  flux  density,  B is  defined  in  the  E frame  as 


1 

1 . , , 1 j ^ 1 

g^  sin  A - hj^  sin  6 1 

i / 

1 / 

3 

f 

= 

Re 

R 

g^  cos  §_+  g'i  sin  6 

/ ^ 

3 

Re 

gj^°  sin  6 + gj^^  cos  6 cos  X + hj^^  cos  B sin  A 


where  R = radius  of  earth  = 2.09029  x 10  ft. 

® 7 

R = orbit  radius  8.  71627  x 10  ft. 

0 = latitude 

A = longitude 


48 


where  if  ^ is  in  webers/m^ 
= -3,04012  X 10"^ 
= -2.1638  X lO"^ 

= 5.7782  X 10”^ 


Latitude  and  longitude  nay  be  calculated  from  the  following  expressions, 


9=0)  cos  ij),, 
o H 


A = 0) 


sin 

^ — - 0) 

o cos  e 


where  oi  = orbital  rate  = 1.454441  x 10  ^ rad/sec  = ^ 

0)  = earth  spin  rate  = 7.272205  x 10  ^ rad/sec 

e 

'('jj  = heading  angle 


The  rate  of  change  of  the  heading  angle  is  given  by 

<l'u  = sin  tan  6 w 
tt  ti  o 


Initial  conditions  may  be  defined  as 

6 * o 

o 

X » o 

o 

= (90°  - i)  , i = inclination  angle  - 63° 

The  flux  density  must  be  transformed  into  body  coordinates  or 
" ^LV/  ^E/LV 

D 
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where 


TE/LV  = transformation  from  E frame  to  local  vertical  frame 


sin  >1;^  - cos  o 

cos  - sin  i|)jj  o 

o o 1 


transformation  from  local  vertical  frame  to  body  frame 
(already  computed  in  program) 


The  resulting  torque  expression  then  is 

!b  = !!b 

Prior  to  launch  the  spacecraft  dipole  moment  will  be  calculated  and 

2 

reduced  to  less  than  .5  anpere-neter  (500  pole-cm)  per  axis.  Thus  if  ^ 
is  given  ampere-meter^  and  is  in  Weber/m^  , is  in  newton  meter.  The 
conversion  factor  between  newton  meters  and  ft.  lbs.  is  .737757  ft-lb/nt-m, 
which  will  be  applied  to  Tg , or 

Tg  (ft.  lbs.)  = . 7375  7 Tg  (nt  m) 

Operating  Range 

Although  it  is  planned  to  measure  the  spacecraft  dipole  moment  prior  to 
launch  and  compensate  to  obtain  values  less  than  500  pole-cm  per  axis. 

Table  F-1,  which  was  taken  from  NASA  SP-8018  indicates  there  may  f.e  a large 
disparity  between  prelaunch  measurements  and  on-orbit  values.  In  magnitude 
alone,  the  values  given  in  I’able  F-1  range  from  a factor  of  .25  to  12.  for 
spin  axis. 
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Table  F-1.  Change  in  Spacecraft  Dipole  Moment 


Spacecraft 

Spin-axis  dipole  moment 
measured  prior  to  trans- 
portation to  launch  site 

Spin-axis  dipole  moment  in 
orbit  (computed  from  motion 
of  spin  axis) 

A-m^ 

(pnle-cm) 

mBSM 

(pole-cm) 

Tiros  II 

+0.1 

(+100) 

+0.9 

Tiros  III 

+0.34 

(+340) 

-0.45 

(-450) 

Tiros  V 

-0.40 

-0.56 

(-560) 

Tiros  IX 

+0.08 

(+  80) 

-0.02 

(-  20) 

ESSA  II 

+0.01 

(+  10) 

+0.10 

(+100) 

ESSA  III 

0.00 

+0.10 

(+100) 

ESSA  IV 

+0.03 

(+  30) 

+0.35 

(+350) 

ESSA  V 

(a) 

(a) 

+0.05 

(+  50) 

ESSA  VI 

+0.17 

(+170) 

-0.15 

(-150) 

NOTE:  Specified  spin  axis  dipole  moment  for  all  spacecraft 
listed  is  0.1  A-m^  (100  pole-cm). 

(a)  Information  is  not  available. 


Table  F-2.  Factors  for  Estimating  Spacecraft  Dipole  Moment  (M) 


Category  of  magnetic 
properties  control 

Estimate  of  dipole  moment 
per  unit  mass  for  non- 
spinninq  spacecraft 

Estimate  of  dipole  moment 
per  unit  mass  for  spinning 
spacecraft 

A-m2/kq 

(pole-cm/lb) 

A-m^/kq 

(pole-cm/lb) 

Class  I 

1 . 'xlO"^ 

(0.45) 

0.4x10'^ 

(0.18) 

Class  II 

3.  xlO"^ 

(1.6) 

1 .4x10"^ 

(0.63) 

Class  III 

10x10“^ 

(4.5) 

4.0x10"^ 

(1.8) 

ah'*  higher 

and  higher 
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Table  F-2,  also  taken  from  NASA  SP-8018  suggests  using  a value  of  1.6  pole 
cm.  per  lb.  of  spacecraft  which  would  suggest  a value  of  1600  pole-cm.  for 
NTS-2.  The  Class  I category  is  for  cases  where  magnetic  torques  are 
comparable  to  other  torques. 

Although  Table  F-1  does  indicate  factor  of  10  variations,  it  appears  that 
these  occurred  for  systems  who  were  tuned  to  finer  than  500  pole  cm.  An 
operating  range  of  + 5000  pole  cm.  is  clearly  unreasonable.  However, 
based  on  Table  F-2,  an  operating  range  of  ± 2000  pole-cm.  would  be  real- 
istic if  slightly  conservative  value. 


I 
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APPENDIX  G 


SOLAR  TORQUE  MODEL 


The  relative  position  of  the  satellite's  surfaces  are  given  in  Figure  0-1 
below,  ^ A 


Figure  G-1.  Relative  Position  of  Satellite's  Surface 


The  solar  force  on  a flat  surface  is  given  by: 


F = -F^N  + F^T 


where : 


F^,  = [(l+aP)cos  ii  + ^ (1-0)  ] P A cos  il- 

N i 

F^  = (1-op)  cos  i'  sin  P A 

- 7 2 

P = Solar  Pressure  Constant  = 1.  x 10  Ib/ft 


(G-1) 
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A = Area  of  the  flat  rsurface 
p = Total  reflected  of  the  incident  light 

0 = That  part  of  p that  is  specular  reflection  (there  is  also  the 
diffused  reflection  6 and  p = (a  + 5 ) 

^ * Normal  unit  vector  to  the  flat  surface 

t = Shear  (tangent)  unit  vector  on  the  flat  surface  given  by: 


(N.l)  N-S 

? = Unit  vector  on  the  sun  LOS, 

ip  = cos  ^ (^,?) 


(See  Figure G-2) 


Figure  G-2.  Unit  Vector  on  Flat  Surface 

The  torque  due  to  (1)  is: 

f “ i xf  where:  (G-2) 

^ » Vector  of  the  center  of  pressure  of  the  surface  (usually  the 
geometrical  center). 
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(G-3) 


The  force  and  the  torque  on  i full  cylindrical  surface  are  given  by: 

I = -F„il  + where: 

N L 

T = Unit  vector  along  cylinder  axis 
^ = Normal  unit  vector  in  middle  of  the  cylindrical  surface  given  by: 

^ T 

t ~ (5.T)  f 

^ = Unit  vector  on  LOS  to  the  sun 

Fj^  = [(1  + od/3)  cos  j + p (1-a)  -1]  P A cos  i|; 

F^  = (1-ap)  cos  ij'  sin  P A 

A = Projected  area  of  the  cylinder  = 2 R H 

R = Cylinder  radius 
H = Cylinder  Height 

i/  = cos  ^ (S.N) 

-7  2 

P = Solar  constant  (1.  x 10  Ib/ft  ) 

The  torque  due  to  G-3  is : 

? = D X ? where:  (G-4) 

S = 2 + rS 

C = Vector  to  the  geometrical  center  of  the  cylinder  (see  Figure  G-3) 
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Figure  G-3.  Vector  to  Geometrical  Center  of  Cylinder 


1 
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TIP  MASSES 


I 


The  total  solar  force  on  a sphere  goes  through  its  geometrical 

center.  Due  to  the  symmetry  each  quarter  (of  the  half  surface  of  the 

sphere  hit  by  the  sunlight)  will  produce  one  quarter  of  this  force  along 
* 

LOS. 


Let  us  divide  this  quarter  into  k equal  slices  by  means  of  k meridians 
angle  Intervals  (Figure  G-A) . 


at 


A $ = 


2k 


Then  let  us  divide  each  slice  into  n spherical  trapezia  by  means  of  n 
parallel  circles  at  Af3  = ^ angle  intervals,  from  the  equator  to  the 
north  pole. 


As  (k  oo)  and  (n  ^ »)  the  spherical  trapezia  can  be  approximated  by 
flat  trapezia. 

The  dimensions  of  the  1th  trapezium  on  a slice  are:  The  two  equal  sides: 

, r = radius  of  sphere  (G-5) 

the  unequal  sides  of  the  trapezium  (its  bases): 

= r^Aij),  r.  = radius  of  ith  parallel  circle 

>^1 
or 


*Other  forces  (tangential)  cancel  out  due  to  the  symmetry  of  the  sphere  and  the 
assunptlon  that  sun  rays  are  parallel  at 


57 


p 


i=l,  2,  where  r = r. 

’ o 

CG  Location 

We  have:  W = R + U (See  Figure  G-5) 

R = (rcos  $ , rsin  4> , 0) 


Figure  G-5.  CG  Location 


(G-6) 

(G-7) 


(G-8) 

(G-9) 
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Figure  G-7 


Figure 


See  also  G-6  and  G-7 


and : 


= r sinC  - r sin  ((1-1)A0) 

= r [sin(lA0)  - sin((i-l)A0) ] 


Substituting  G-11  and  G-13  into  G-12  we  obtain  the  ^ vector  of  the 
1th  trapezium 


1) 


and  are  given  by  G-14,  G-15,  and  G-16 

2)  X is  given  by  (5a)  and 

3)  $ = (j-l)Ai^  + , j = 1,  2,  ....  (k-1) 

= (j-0.5)A^ 

j = jth  slice  (j  = 1,  2,  . . . , k) 

i = ith  trapezium  on  this  slice  from  equator  to  north  pole. 
(i=l,  2,  ...,n) 


C.  The  normal  unit  vector  to  the  ith  trapezium  is: 


-> 

= 

CG 

1 

CG 

where 

-¥ 

CG  is  given  by  G-17  and  the  norm: 

CG 

= r or 

CG 

=~\J  [xcos 

4’]^  + [xsin  + [c^  + rsin(  (i-l)A6) 

+ [c  + rsin((i-l)Ae) ] ^ 
z 


(G-18) 
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I 

GRAVITY  GRADIENT  RODS 

Let  S be  the  length  of  the  rod.  The  rod  is  a full  cylindrical  surface 
and  when  it  is  curved  due  to  thermal  bending  we  can  approximate  it  by 
N equal  small  cylinders,  each  having  height  S/N  and  which  we  assume  to 
be  rlg^it  cylinders.  See  Figure  G-9. 


62 


Figure  G-9 

Let: 

X Y Z - Body  Frame 

X y z - Rod  Frame  (same  as  XYZ  at  nominal  position) 

R = Radius  of  curvature  of  rod. 

We  have:  a = ^ , N 

,i  = 0,l,2,  (N-1) 

0 < 't'  < a 

From  Figure  (;-9  we  dec. '.'mine  the  normal  unit  vector  to  the  ith  sub 
division  in  Rod  Frame. 


= .(-Rcos  4> , 0 , -Rsin  4>) /R  = (-cos  0,  0,  -sin  <I>) 
where  $ is  given  by  G-19. 


Let  be  the  transformation  from  Rod  Frame  to  Body  Frame. 
RB 

Then  N = T N is  the  normal  unit  vector  in  Body  Frame. 

RB  X 


X 

X 


We  must  determine  T , i.e.,  the  three  Euler  angles. 

Kd 

Let  S = (Sj^,  S2 , S^)  be  the  sun  unit  vector,  and  its  projection 
vector  on  the  XY  plane.  Then(  see  Figure  G-IO)  we  have 


(G-20) 


and  we  have  the  Euler  angles  for 
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Figure  G-10 
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Let  us  deterraLne  the  unit  normal  vector  N of  the  ith  rod  element  in 
Body  Frame  from  (corresponding  one  in  Rod  Frame)  and  the  known  Euler 
angles  (6,  li/,  9)  given  by  G-20. 


We  had:  = (-cos  4>,  0,  +sin  it),  $ = (i-.5)A9 

or  N,  (x  , 0,  z ) i = 1,  2,  N 

l o o »»» 

= (x^,  0,  z^)  in  xyz  rod  frame 


(G-23) 


i)  Find  in  x'y ' z ' f C call  it  N2)  . See  also  Figure  G-11. 
For  this  project  x^  on  x'-y'-z'  axes  actually  x'-z'-axes 
(since  the  y'  projection  is  always  zero  because  is 
always  on  the  x'z'  plane). 


This  gives:  (x^  cos  a,  0,  ±x^  sin  a) 

Project  z on  the  x'z'-plane:  (+z  sin  a,  0,  z cos  a) 


= (x  cos  a + z sin  a,  0,  ± x sin  a + z cos  a) 
2 o 0 0 o 


(G-24; 


ii)  Find  ^2  XYZ  Body  Frame  (call  it  iJ) 

For  this  we  need  only  project  (x  cos  a + z sin  a)  on  the 

o o 

X,  Y axes  since  Z and  z'  are  identical. 


^ = ((x 

cos 

a + 

z 

sin 

a)  cos  6,  (x  cos  a + z sin  a)  sin  9, 

(G 

0 

0 

0 0 

± X 

sin 

a + 

z 

COS 

a) 

0 

0 

where : 


x^  = -cos  $=(i-0.5)Acp,  i=l,  2,...,N 

z =±sin$;Acp=-^ 
o ^ N 
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where:  the  top  sign  applies  to  the  rod  on  -Z  axis  and 

the  bottom  si^  applies  to  the  rod  on  +Z  axis. 


Finally  we  must  determine  the  CENTER  vector  in  body  frame  at 
the  ith  element  of  the  rod  as  well  as  its  AREA. 


From  Figure  G-12  we  have 


CENTER  = R (S  + I)  + R S 
curv 


Figure  G-12 


O 1 u ^ 

where  f = ( , — ,0),  ^ = (S  , S_,  0.)  unit 

ISII  Hll  ‘ 


R = Rod  Radius 

^ is  given  by  G-23  and  R is  the  input  radius  of 
^ curv  ^ 

curvature  of  the  rod  due  to  the  thermal  bending. 
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Definition 


Definition 


AREA.  = Acp  * R * WR0D 
1 curv 

WR0D  = Rod  Diameter 

When  R = i.e.,  straight  rod  the  formulae  above  must 

curv 

be  modified  as  follows: 
e = tan  ^ (82/5^) 

a = 0 
Acp  = 0 
$ = 0 

= (-1,  0,  0) 

S = (-cos  6 , -sin  6 , 0) 

► 

CENTER  = (0,  0,  +SR0D/2  + RBIAS) 

SR0D  = S (Rod  Height) 

AREA  = SR0D  * WR0D 


SOLAR  PANELS 

0 = angle  between  the  panel  normal  and  LOS  to  the  sun 

P 

and  is  positive  if  rotation  is  CCW  from  LOS  to  the  sun, 
when  standing  at  the  origin  of  the  body  frame  and  look- 
ing at  +Y  axis. 

Panels  are  at  nominal  position  when  their  normal  unit 
vector  is  [5  = (1,  0,  0) 
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* • — * Nominal  position 

Figure  G-13 

-1 

From  Figure  G-13  we  have:  <)>  = cos  ( 

IIS  II 

S = (S^,  S2,  S^)  unit  vector  along  LOS. 

Then  i()  = <}>  + 9 and 


“V 

N » (cos  ij),  0,  sin  ip)  is  the  normal  unit  vector  to  the  panel 


Figure  G-14 


From  Figure  G-IA  we  see  that: 

CENTER  = CP  + D2  + Ar 
where : 

CP  is  input 

Ap  = Solar  pressure  shift  vector  from  the  panel  geometrical  center 
(given  in  Body  Frame  when  panel  is  at  nominal  position)  and  is 
input. 

2 = (0,  d/2.  0) 
d = length  of  panel 
h = width  of  panel 
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NTS-2  GRAVITY  GRADIENT  ROD  THERMAL  ANALYSIS 
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I 
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The  BI-STEM  boom  configuration  and  2-D  thermal  model  are  shown  in  Figure  H-1 . 

The  NRL  specification  quotes  a maximum  solar  absorptivity  of  0.15;  no 

emissivity  data  w.is  available  - a nominal  value  of  0.05  was  assumed. 

Preliminary  calculations  showed  that  compared  with  peripheral  conduction. 

internal  radiation  could  be  neglected.  Radiation  was  Included  between  the 

overlapping  sections  of  the  boom  elements.  A hot  case  solar  flux  of 
2 

455  Btu/hr.  ft.  was  assumed.  Earth-albedo  and  IR  fluxes  were  excluded, 
providing  a major  analysis  simplification.  Inclusion  of  these  terms  would 
result  in  a slight  reduction  of  the  predicted  temperature  gradients. 

Note  that  the  model  is  two  dimensional,  i.e.,  axial  temperature  gradients  in 
the  boom  are  not  considered. 

Temperature  Predictions 

Figures  H-2  to  H-7  show  steady  state  temperature  for  two  solar  orientations  ai.r 
for  various  values  of  a,  e.  Figure  H-2  is  for  nominal  a,  e and  indicates  a 
gradient  of  only  2°F  from  sunside  to  backside.  Figure  H-3  shows  the  effect  of 
severe  degradation  in  both  a and  e:  the  gradient  increases  to  4.l'^F.  (Note 
that  the  mean  boom  temperature  is  unchanged  becaijse  the  a/e  ratio  is  the  same.) 
Figure  H-4  shows  the  effect  of  increasing  the  emissivity  only;  although  the 
mean  boom  temperature  drops  significantly,  the  gradient  is  still  only  2°F . 

Figures  H-5,  H-6  and  H-7  stio’’  temperatures  for  the  same  property  comblnatic ns  but 
for  a different  solar  orientation.  Hot  side  to  cold  side  gradients  are 
unchanged. 
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DIAMETr-H  0.5" 

ELEMENT  [HIOKNESS  1.002" 
MATERIAL  RcCu  Alloy  No.  ^5 


Figure  H-1.  Bi-Stem  Boom  Configuration  and  Thermal  Model 


f 
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figure  H-8  shows  the  transient  response  for  entering  and  leaving  a 55  minute 
period  in  the  earth's  shadow.  Initial  temperatures  were  assumed  to  be  those 
shown  in  Figure  H-5;  at  time  zero  the  solar  flux  was  removed  and  the  tempera- 
ture histories  computed  for  55  minutes.  At  this  time  the  sun  was  switched  on 
again  and  the  response  confuted  for  a further  200  seconds.  The  results  in 
Figure  H-8  show  that  the  2°F  gradient  decays  within  10  seconds  of  eclipse  and 
re-establishes  within  10  seconds  of  re-entering  sunlight.  In  actuality,  the 
change  would  be  slower,  since  the  transition  from  sunlight  to  darkness  (and 
vice  versa)  would  take  a minute  or  so  rather  than  being  a step  function  as 
assumed  in  the  analysis. 

Boom  Deflections 


To  get  an  estimate  of  tip  deflection,  consider  the  relative  expansion  of 
"strips"  on  the  sunside  and  shadeside: 


ft  = 1 _ ^ + gg- 
® R " R + d 

l.e. 

R + d _ + 6 i 

R 


d X g- 
S I 


BO 


But 

SZ  = Bx  at  xH 

where  B is  the  coefficient  of  expansion 

* R 1_ 

• • ^ B X AT 

For  Beryllium  Copper  Alloy  No.  25, 

B = 9.3  X 10”^  ln/in°F 

Thus  for  a AT  of  2°F,  and  d = 0.5" 

R = = 26,881.7  inches  = 2240  ft. 

9.3  X 10  X 2 

The  tip  deflection, 

6 = R (1-cos  6) 

' 

« t 6 = 0.83  ft.  (For  2 F gradient) 

[For  comparison,  the  formula  on  P.13  of  NASA  CR-2516  gives  a static  tip  I 

deflection  of  1.05  ft.  for  this  geometry.] 

i 

It  had  been  planned  to  generate  a STARDYNE  model  of  a finite  length  of  the 
boom  and  impose  the  computer  cross  section  temperature  profile  and  determine 
displacements.  In  view  of  the  small  gradients  and  tip  deflection,  this  is 
not  presently  considered  necessary. 


Jo  jbfiivajdwji 
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Figure 


Implementation 

The  thermal  model  was  set  up  to  run  a series  of  cases  for  solar  incidence 
angles  verying  from  0°  to  80°.  The  configuration  was  simulated  with  "nominal" 
thermal  properties  Temperature  gradients  are  summarized  in  the  table  below 
and  plotted  in  Figure  H-9.  Full  thermal  maps  are  given  in  Figures  H-10  through 
H-18. 


Table  H-1 . Temperature  Gradients 


SOLAR 

INCIDENCE 

0 

10° 

0 

o 

CM 

30° 

o 

0 

50° 

60° 

O 

0 

o 

o 

00 

S unside 

Temp.  °F 

251.59 

248.42 

240.12 

225.89 

205.06 

176.40 

137.52 

83.26 

-1.51 

Shades  ide 

™ 0-^ 
Temp . F 

249.53 

246 . 39 

238.18 

1 

224.11 

1 

203.49 

175.07 

136.49 

1 

82.56 

-1.86 

T °F 

2.06 

2.03 

1.94 

1.  78 

1.57 

1.33 

1.03 

0.70 

0.35 

Static  Tip 

Deflection 

(Feet) 

1.0 

0.985 

0.942 

0.864 

0.762 

0.646 

0.50 

0.34 

0.17 

The  boom  static  tip  deflections  shown  assume  that  a 2°F  gradient  gives 
approximately  a 1 ft.  deflection. 
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▼ 
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Figure  H-12 
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50°  INCIDENCE  ANGLE 

I SUN 


176.4 


Figure  H-15 
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Summary  and  Conclusions 

For  nominal  external  thermal  properties  (a  = 0,15,  e = 0.05)  the  steady  state 
temperature  gradient  across  the  boom  is  predicted  to  be  only  2°F.  This 
would  produce  a deflection  of  about  1 foot  at  the  boom  tips.  With  degraded 
thermal  properties  (a  = 0.3,  e = 0.1)  this  gradient  would  increase  to  4°F 
and  the  tip  deflection  to  2 feet. 

Transient  thermal  response  was  evaluated  for  step  function  sun  on/sun  off 
conditions.  The  results  show  that  the  2°F  sun  side  to  back  side  gradient 
is  established /or  decays  very  rapidly  (within  10  seconds) . 

This  analysis  suggests  that  deflection  of  the  booms  as  a result  of  thermal 
gradients  should  not  be  a problem.  The  effect  on  the  spacecraft  control 
system  could  be  assessed  by  imposing  a 1 - 2 ft.  step  displacement  on  the 
boom  tips.  If  problems  occur  the  displacement  could  be  ramped  over  time 
periods  of  10  seconds  to  1 minute  instead  of  assuming  a step  function. 
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appendix  I 


GRAVITY  GRADIENT  ROD  FLEXIBILITY  MODEL 

I-l  NTS-2  GRAVITY  GRADIENT  ROD  DYNAMICS 

This  appendix  contains  details  of  Stardyne  simulations  of  a rigid  satellite 
with  two  flexible  gravity  gradient  rods.  These  simulations  showed  that 
modeling  each  rod  with  one  flexible  beam  and  one  tip  mass  gives  questionable 
results.  Using  two-beam  elements  or  four-beam  elements  gives  almost 
identical  results.  Therefore,  a two-beam  model  was  further  explored  for 
addition  to  the  orbit  simulation  model  (GPS) . (See  Section  1-2) . 

Contents 

• Summary  Table  - Key  Results  (Table  I-l) 

• Deformed  Geometry  Sketches,  for  modes  responsive  to  solar 
array  drive  torques  (Figures  I-l,  1-2,  and  1-3) 

• Math  Models: 

2 Nodes /Rod  (1  beam/rod),  original  (Figure  1-4) 

2 Nodes/Rod  (1  beam/rod) , updated  (Figure  1-5) 

3 Nodes /Rod  (2  beam/rod)  (Figure  1-6) 

5 Nodes /Rod  (4  beam/rod)  (Figure  1-7) 

• Modal  Data  Summaries,  Tables  i-6  and  1-7 

• Mass  Property  Calc,  for  Dynamic  Rod  Models 


TABLE  I-l.  SWMARY  TABLE— RESULTS  OF  FLEXIBLE  GRAVITY  GRADIENT  ROD  STUDIE 


FOR  DYNAMIC  MODEL  WITH  2 NODES /ROD 


11 

I Q ' 
' T 


NODE  6 y 


NODE  2 


.45  lbs.  0 center  of  upper  mass 

Rod  Weight  = .90  lbs . 

.45  lbs. 


19.75" 

ADDITIONAL  INERTIA  FOR  NODE  31: 

AIxx  = (-^5)  2 flj.86^  + 11.14^)j 

518.5 

= 466.67  Ibs-in^ 

AIyy  = (-45)  (2)  (28.34^+  22.95^) 

= 1196.87  Ibs-in^ 

41^2  = (.45)  (2)  [19.75]^ 

= 351.06  Ibs-in^ 

AW  = .90  lbs.  total  added  to  Node  31 
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updated  dynamic  model.  2 NODi-S/ROD  (TIP  and  SPACECRAFT  END) 


For  cantilever  beaa  of  mass  and 
tip  mass  m 

Equivalent  tip  mass  weight  = m + 

Ref.  MyKlestad,  N.O,,  "Fundamentals 
of  Vibration  Analysis,"  McGraw-Hill, 
1956,  p.  40. 

Above  is  based  on  Raylei^-Ritz 
energy  method. 

Therefore,  transform  % beam  weight  to 
end  of  beam  but  do  not  transform  h beam 
inertia,  since  this  is  accounted  for  in 
the  H m^  approximation. 

Node  11 

AW  = .225  lbs. 

AIxx  = lyy  ° ^-8125)^  = 

3.2  7 Ibs-in^ 

Node  21 

AW  = .225  lbs. 

AI^  = AI^y  = .225(2)^  = 0.9  Ibs-in^ 


I 

I 

I 

I 

I 
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N ode  31  (Spacecraft  c.g.) 


Assume  3/4  of  rods  attached  and  full  inertia  transformed  to  base. 
AW  = 2 X .675  = 1.35  lbs. 

From  upper  rod 

Alyy  = 1/3  W + W(19.96^  + 19.75^) 

= l/3(. 675)540^  + .675(  ) = 

= 65,610  + 532.2 

2 

= 66,142.2  lbs- in  ^ from  upper  rod  only 

From  lower  rod 

AIyy  = 1/3  W 8.^  + W(11.14^  + 19.75^) 

= 65,610.  + .675  ( ) 

= 65  ,610.  + 347. 

= 65,957.  Ibs-in^.^ 

Total  Iyy  added  to  Node  31 
= 66,142.2  + 65,957 

‘'yy  total  ■ 

, Tor 

From  upper  rod 

AIyy  = 1/3  W + W 19.96^ 

= 65,610  + .675  (19.96)^  = 

I u, — _ j 

I 268.9 

, = 65,878.9 
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From  lower  rod 


AIyy  = 1/3  W + W(11.14)^ 

= 65,610.  +,  .675  (11.14)^  i 
83.  76 

= 65,693.76 

"xx  TOTAL  - 131.572-’  H>A-Ln" -• 

41^2  = 2 X W . (19. 75)^ 

= 2 X .6  75  X 19.75^ 

AI22  = 526.584  Ibs-in^  
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3 NODES  PER  ROD 


Program  will  automatically  lump  rod 
masses  to  Nodes 


2,  4,  6 (upper  rod) 

12,  14,  16  (Icwer  rod) 

Will  manually  transform  weights  at 
dependent  nodes  2 , 6 , 12  and  16  to 
adjacent  Independent  nodes  of  rigid 
systems . 

Node  11; 

AW,,  = .225  lbs.  (W  , W , W ) 

11  x’  y z 

AI  = AT  = .225  (3.8125)^ 

A I 

=3.27  Ibs-in 


Node  21: 

AW_,  = .225  lbs.  (W  , W , W ) 

21  x’  y z 

AI^  = AI^  = .225  (2.0)^  = .90  Ibs-in 


Node  31: 


AW^^  = 2 X (.225)  = .450  lbs. 

(W  , W , W ) 

X ’ y ’ z 

AI^  = .225  (19.96)^  + .225  (11.14)^ 

= 117.56  Ibs-in^ 

Aly  = .225  (19.96^+  19.75^)  upper 

rod 

+ .225  (11.14^  + 19.75^)  lower 

rod 

= 293.09  Ibs-in^ 


A I„  = .225  X 2 X (19.75)' 

^ 2 
= 175.53  Ibs-in 


3 NODES /ROD 


MOVE  Z MASSES  TO  SPACECRAFT 


1 


f 


Program  generates  rod  wei^ts  at 
nodes  2,4,6 
12,  14,  16 

Must  manually  transform  weight  to 
independent  nodes  of  any  rigid  systems. 

Node  11:  (from  Node  6) 

AW,,  = .225  lbs.  W , W 
11  X ’ y 

N ode  2 1 : (from  Node  16) 

AW„,  = .225  lbs.  W . W 
21  X z 

Total  rod  + tip  mass  weight  in  Z 
direction  = 7.5  + .90  = 8.4  lbs. 
Transform  to  Nodes  2 and  12,  which  are 
part  of  rigid  systems.  Hang  Z weight 
on  new  Nodes  at  ends  of  stiff  elastic 
beams  tied  to  Node  31.  Node  23  at 
same  location  as  Node  2,  Node  24  at 
same  location  as  Node  12. 


For  beam  from  Node  23  - 31 


Want  f 


cantilever  beam 
50.  H 


n 


,,  2 _ ...2  K 386.4 

”n  - 50) 


16.8  X (IOOit)' 
386.4 


= 4291. 


4291  (19.75)  ^ 

3 X 18.5  X 10^ 


Thus  total  c.g.  is  correctly  modeled 
and  inertia  of  Z weight  contribution 
is  retained. 


E = 18. 5 X 10 


b = 12kl 


0.030  1 6 for  b = d b = . 775 

I = 1/12  b^ 
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Figure  1-7 
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ADDITIONAL 

INERTIA  FOR 

nod;.  31: 

2 

O 

AI 

= .1125  (2) 

19. > 6^  + 

11. 14*- 

XX 

9 

= 116.67  Ib- 

, L 

in 

AI 

= .1125  (2) 

2 

28,  + 

22.95^ 

yy 

2 

= 299.22  Ib- 

in 

AI 

zz 

= .1125  (2) 

19.75  ^ 

= 87.764 

= .1125  X 2 = .225  lbs  to  Node  31 


AW 


TABLE  1-2.  STARDYNE  MODAL  DATA  1 BEAM  PER  GRAV.  GPJiiD.  ROD.  ROD  El  = 3000.  LBS-IN 
SOLAR  ARR^\YS  ASSUMED  RIGID  BUT  FREE  TO  ROTATE  ABOUT  X2  AXIS. 


ape 


TABLE  1-3.  STARDYNE  MODAL  DATA  1 BEAM  PER  GRAY.  GRAD.  ROD.  ROD  El  = 3000.  LBS-IN 
SOLAR  ARRAYS  ASSUMED  RIGID  BUT  FREE  TO  ROTATE  ABOUT  X2  AXIS  (updated)** 


(H')K  riir.  I'SNSr  \\  modes  iin  tapo  X7b5) 


TABLE  1-4.  STARDYNE  MODAL  DATA  2 BEAMS  PER  GRAV.  GRAD.  ROD.  ROD  El  = 3000.  LBS-IN“ 
SOLAR  ARRAYS  ASSUMED  RIGID  BUT  FREE  TO  ROTATE  ABOUT  X2  AXIS. 


2 


HQR  run  TSNSTLi,  nodes  oi  t.ipt*  X417ft 


TABLE  1-5.  STAFDYNE  MODAL  DATA,  4,  BExWS  TER  CRAY.  GRAD.  ROD.  ROD  El  = 3000.  LBS-IN 
SOLAR  ARRAYS  ASSL-MED  RIGID  BUT  FREE  TO  ROTATE  ABOUT  X2  AXIS. 


NOTES:  *Spact.*cr.ift  c.g. 

fRun  No.  TSNSTA5,  Moes  Tapt-  X3381) 


TABLE  I-’ 


SPACECRAFT  PARAMETER  VALUES 


I 

J 


Mass  Properties 

(orbi tal  configuration) 

Total  weight  (lbs)  = 910 

Orbit  correction  - 75t  hydrazine  usage 

Spacecraft 

moments  and  products  of 

i ne  r t i a 

, solar  arrays 

deployed. 

gravity  gradient  rods  not 

deployed  (slug  - ft  ) 

^\.x  ° 

208.  89 

^XX  " 

0 

lyy  . 

167.46 

^XZ  " 

0 

^ZZ  “ 

78.98 

4z  = 

-3.04  = 14082.73  Ibs-ln^ 

spacecraft 

moments  and  products  of 

inertia 

with  arrays  and  rods 

deployed  (slug  - ft*") 

^XX  " 

209  = 968188.3 

Ibs-in" 

0 

lyy  - 

1998  - 9.2557  x 

10^ 

'xx-  = 

0 

^ZZ  " 

1909  = 8.8434  x 

10° 

^YZ  ' 

-3.04 

Solar  arrav  moment  of  inertia  about 

arrav 

axis  of  rotation 

“ 3.28  slug  - ft^  (total  1 

for  two  arrays)  = 

15194.5  lbs-in‘. 

Spacecraft 

weights,  center  of  gravity,  moments  and  products 

of  inertia 

Cslug-ft*")  and 

liy  drazine 

usage 

for  orbit  connection 

in  stowed 

condition. 

Total 

C (in) 

‘zz  ^YZ 

Hydrazi ne 

WT  (lbs) 

^XX 

Usage 

9U9 

25.75 

124.  12 

73.91 

88.78  -3.31 

0 

936 

25.58 

122.55 

71.92 

88.  34  -3.31 

25- 

923 

25.44 

120.94 

70.08 

88.05  -3.31 

50 : 

910 

25.  33 

119 . 34 

68.33 

87.85  -3.31 

752 

89  7 

25.26 

117.75 

66.69 

87.  76  -3.  31 

lOOX 

1 

I 


TABLE  1-7.  INPUT  TO  RUN  G2  MASS  INERTIA 
PROPERTIES  COMPUTER  PROGRAM 


NTS-2  with  rigid  deployed  si'lar  arrays  less  elastic  grav.  grad,  rods 
Total  properties  from  G.E.  report,  NRL  Table  4-1.  c.g.  from  top  of  vehicle. 


Total  Weight 

2 

910. 

1 Start  with 

Total  c.g. 

19.86 

> total  rigid 

Total  Mi. 

9255695.04 

8843404.32 

968188. 32 

' spacecraft 

Total  Prod. 

00 

O 

1 

Uprmass  Weight 

2 

-7.5 

Uprmass  c.g. 

19.75 

-723.8125 

Uprmass  Mi. 

-48. 

-48. 

-48. 

Upr  rod  Weight 

2 

-.90 

Subtract  elastic 

Upr  rod  c.g. 

19. 75 

-360. 

rods  and  tip 

Up  r rod  Mi . 

- 38880.03 

-38880.03 

-.0558 

masses 

Ijwrmass  Weight 

2 

-7.  5 

) 

Ijwrmass  c.g. 

19.75 

75  3.0 

Lwrmass  Mi . 

-27. 

-27. 

-27. 

Lwr  rod  Weiglit 

2 

-.90 

Lwr  rod  c.g. 

19.75 

391 .0 

Lw r rod  Mi  . 

-38880.03 

-38880.03 

-.0558 

J 

Endgroup 

E!ndcase 

Stop 

$ END.IOB 

S*$DIS 

$*$DIS 


RESn.TS  for  Node  31,  Flex  Rod  Models 

Wt.  - 893.2 

c.g.  * 

X **  - . 3714  in  . 
y - 0.0 
z = 19.957  in. 

Inertl.vi  .ibout  c.g. 

^XX  ’ lhs-in‘ 

lyy  - 325980  - 15,195 
''  2 

■ 310,785  lbs. -in  for  solar  arrays  free  (.  turn  .il-.- .i  ' ive  axis 
^ZZ  “ 961440.  lbs-in‘ 


TABLE  1-8.  INPROP  RUN  G2  MASS 


COMPUTED  VALUES 

NTS-2  with  rigid  deployed  solar  arrays  less  elastic  grav.  grad  rods 


RESULTS  FOR  THE 

TOTAL  SYSTEM 

Wt.  “ 

0 . 89  320E 

03  +- 

0. 

Pounds 

c. 

.g.  Loc. 

-0.  37147F: 

00  +- 

0. 

Inches 

0. 

0. 

Indies 

0.199 5 7E 

02  +- 

0. 

1 nches 

-0.33180E 

03  +- 

0. 

Pound  - 

Inches 

0. 

0.  ” 

Po'jnd  - 

Inches 

0.17826E 

05  +- 

0. 

Poind  - 

Inches 

Square 

> 

0.11007E 

07  +- 

0. 

Pound  - 

Inches 

square 

Square 

- 

0.68185E 

06  +- 

0. 

Pound  - 

Inches 

square 

Square 

- 

0.96156E 

06  +- 

0, 

Poind  - 

Inches 

square 

Square 

- 

-0. 14083E 

05  +- 

0. 

Pound  “ 

Inches 

sq  uare 

Square 

- 

0. 

+- 

0. 

Poind  - 

Inches 

square 

Square 

- 

-0.48744E 

04  +- 

0. 

Pound  - 

Inches 

square 

Mass  mom. 

of  inertia 

about 

c.g. 

Square 

^XX 

0.74494E 

06  +- 

0. 

Pound  - 

Inches 

square 

Square 

lyV 

0.  3259 8e 

06  +- 

0. 

Poind  - 

Inches 

square 

Sq  uare 

^ZZ 

0.96144E 

06 

0. 

Pound  - 

Inches 

sq  uare 

Square 

-0. 14083E 

05 

0. 

Pound  - 

Inches 

square 

Square 

^z 

0. 

+- 

0. 

Pound  - 

Inches 

square 

Square 

0.174  74E 

04  +- 

0. 

Pound  - 

Inches 

square 
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1-2  DERIVATION  OF  INFLUENCE  COEFFICIENT  MATRIX,  INCLUDING  CROSS-COUPLING 
EFFECTS  DUE  TO  ROD  CURVATURE 

Gravity  Rod  Simulation  with  Two  Flexible  Beam  Elements 

/ (Linear  Theory) 

I 

Based  on  Stardyne  model  vibration  analysis  of  NTS-2  with  rigid  solar  arrays 
I and  flexible  gravity  gradient  rods,  the  following  flexible  rod  model  with 

two  beam  elements  of  equal  length  was  found  to  give  very  good  results. 


The  influence  coefficient  matrix  for  the  above  case  can  be  calculated  using 
standard  beam  deflection  equations,  shewn  on  the  next  page.  Thus 


{ 
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Under  dynamic  loading  cases,  the  net  applied  force  at  each  panel  point  will 
be 

i external^  1 i 

where 

F could  be  due  to  solar  pressures  or  the  coirponent  of 

external 

gravity  gradient  torque  (forces)  normal  to  the  rod. 

From  Roark,  "Formulas  for  Stress  and  Strain,"  4th  Ed.,  the  following  beam 
equations  are  presented. 


TABLE  1-9.  SHEAR,  MOMENT,  AND  DEFLECTION  FORMULAS  FOR  BEAMS; 

REACTION  FORMULAS  FOR  RIGID  FRAMES 

Notation:  W = load  (lb.);  w = unit  load  (lb.  per  linear  in.).  M is  positive 
when  clockwise;  V is  positive  when  upward;  y is  positi’/e  when  upward. 

Constraining  moments,  applied  couples,  loads,  and  reactions  are  positive  when 

acting  as  shown.  All  forces  are  in  pounds,  allmoments  in  inch-pounds;  all  ^ 

deflections  and  dimensions  in  inches.  6 is  in  radians  and  tan  0 = 6 

I 

i 


I 


STATICALLY  DETERMINATE  CASES 


1. 

•adtrig,  Supju'ft,  and 

Ri  dv. ' Ion-  1 uiiU  R 1 , 

Bending  MoiTiini  M and 

D<.  1 1 ■ 1 ■•n  V , ^‘.ix  ir.u'  .it!-  t i , • , 

Rv  fcren^.*  N iitiher 

rt  K-..1  >bi..  jr  V 

Maximum  He-nding  Momint 

.»n.]  I r.d  bio;-. 

'• 

RJ  ■ +K 

't  * -W.\ 

V 

V - -W 

'IdxM  ■ -W  .il  K 

M.„  • - ! „ A 

; 

f 

. - a-  A 

R.)  « +W 

fA  to  H)'1  - II 

1 Fl'--''  • ■'  ■ 

¥ 

0 

..  if 

(,\  to  BU’  » 0 

(B  to  ' )M  » -K(x-b> 

(B  lo  1 IV-  ' (s-  h,  ■ - .w-lv  -■  . , 

(B  to  C )(  « -W 

■ -Wa  .It  1 

May  y . - ^ 

- ‘ ' !^f  ■' 
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Note  that  no  external  beam  t.irque  loads  are  included. 


Definitions : 

e.  . = deflection  of  point  A due  to  unit  load  at  point  A 

e^  - = deflection  of  point  A due  to  unit  load  at  point  B 

^2  2 ~ deflection  of  point  B due  to  unit  load  at  point  B 

2^'^  reciprocity  theorem. 

^2  I ~ deflection  of  point  B due  to  unit  load  at  point  A 


I 

I 

I 
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By  Case  1 


3 3 

_ 1 1 5. 

“ 3 El  3(1-0) 


Rotation  at  point  b,  (by  Case  1) 


= u Pi: 


Deflection  at  point  A due  to  load  at 
point  B is  then 


y^  = yb  + 6 • £ 


'1,2  6 El 


From  first  calculation  for  e 


Numerically,  for  straight  r.ids,  NTS-2 
r41,A72.  12,96oT| 

[e]  = inches /lb. 

_12,960.  5,184. 

Note : 

Define  total  rod  panel  point  deflection  as  follows 
^TOT  ^elastic  ^rigid  body  ^thermal  ^non-linear 


Define  inertia  loads  in  terms  of  y^  ^ , 

^ total 


V *v  “^y  "^y 

total  -^elastic  rigid  body  ■^thermal  •^non-linear 


y , j is  determined  by  means  of  influence  coefficient  matrix  above, 
elastic 


Similar  expression  will  also  govern  elastic  beam  bending  in  the  other 

beam  bending  axis.  That  is,  for  uniform  El  = El 

“ XX  yy 
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Now  examine  effect  of  sliglit  bending  of  rod  due  to  out-of-straightness, 
thermal  deflections,  or  other  sources,  causing  an  equilibrium  position  of 
the  rod  with  a slight  curvature. 

linear  elastic  cross-axis  coupling 


at  Node  A 

Xq^  = R - R cos  6 = R( 1 “ cos  6) 
at  Node  B 

Xpg  = R(1  - cos  (j)) 

Choose  positive  x in  direction  toward  sun,  resulting  thermal  bending 
is  away  from  sun. 

New,  a force  perpendicular  to  x-z  plane  will  result  in  a linear  y direction 
deflection,  as  before  plus  additional  deflections  due  to  twisting  about 
the  z axis . 


123 


In  3-D  View 

For  torsion 
TL 

® JG 


T = torque 
L = beam  length 
J = torsion  section  stiffness 
G = shear  modulus  of  material 


In  Side  View,  in  X-Z  Plane 
~ ^OB 


base 


OB 

i 


For  loading  into  the  plane  of  paper, 

in  the  y direction,  no  torsion  exists 

in  the  beams  if  d-  = * • 

A B 

In  upper  beam  segment,  no  torsion 

exists  due  to  F , 
yA. 

Pure  bending  about  point  B is  about  an 
axis  to  upper  beam,  line  A'  - A’ . 
This  moment  at  B about  A 'A'  must  be 
resolved  into  a moment  about  axis 
B'-B  and  a torque. 


XoA 


= F • z sin 


for  small  (j) 


, ^A  ~ ’^B  ^OB  „ 

= z sin  F , = (x„,  - 2 Xr,„)F  , 

Z Z yA  OA  OB  yA 
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At  Node  B 


vA 


Twist  of  lower  beam  Is  0 = 

5 

e = [(""OA  ~ ^ "^OEoli 
B JG  yA 


Ti 

JG 


Additional  deflection  in  y direction  at 
Node  B = 0 


However,  twist  6^  about  Icwer  beam  will  cause  some  additional  y deflection 
at  upper  node. 

Arm  = i sin  (4),  - 4>„) 

A D 

■ '«A  - >B> 


= i 


^A  ■ ^OB 


OB 

i 


Deflection  at  end  of  arm  (Node  A)  is  9_  • arm 

D 


JG 


yA 


for  = 36 


OB 

S>.  = 360 . 

JG  = 25 


above  term  is 


(36  - 18)  360. 

25. 


= 4665.6  In/lb. 


Checks  very  well  with  3 node  Stardyne  static  model.  Low  for  5 Node  static. 
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Coupling  with  Acceleratlor.  along  Z Axis  (Rod  Axis) 


As  shown,  a posltl\’e 

force  F will  tend  to 
z 

strai^ten  the  rod. 

A negative  oriented  Z 
force  will  Increase  the 
rod  deflection. 


Buckling  check: 

From:  Myklestad,  N.O. , "Fundamentals  of  Vibration  Analysis,"  1956, 
p.  28. 


cr 


/ T 2 

^ L mg 
0 

-■  El 


for  F = ng  dovn 
z 


P = column  buckling  load 
m = ( 7.5  + .45)  7386.-  Ibs-sec'^ /In . 


g = dynamic  load  factor  acting  dcxgnward  due  to  reactions  from 
solar  array  drive  torques,  inches /sec”. 

2 , 

mg  = - El /4  L“ 


El  g 


=allow 


c -■  (3000.)386.4  2 

- = ^ r = o.ol  in /sec 


4 L" 


(720)' 
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Coupled  X-Z  Beam  Deflectlor, 

Mnear  theory  will  be  used.  Exact  buckli  ig  solution  involving  theory 
of  elastic  stability  is  not  suitable  for  application  at  hand,  especially  since 
coluim  loadings  are  very  low,  by  prev.  calculations. 


; 


sin  0 = 


R 


^OA 
cos  0 


R - R cos  0 = R( 1 - cos  0) 


1 


R 


tan  0 

9 = 


Xqb  = R(1  - cos  -1) 


No  closed  form  solution  available  to  solve  for  R in  terne  of  I, 
Must  Iterate.  ^ 


f 
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Assume 


'^OA’  OB 


are  kncwn,  estimated,  or  previously  calculated. 


Resolving  vertical  force 


F = F 
xB  zA 


sin  ^ 


transferred  to  point  B,  lateral  component 
zA  ‘ ‘*’B 


is 
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2-X  Coiroled  Deflections,  cmitlnued 


Deflections  of  point  B (elastic) 

1 2 
, F „ - M „ 1 

. = I ^ + i 

^ 3 El  2 El 

1 ’’zA-  1 ♦a 

^ “ 3 El  2 El 


Rotation  of  point  B (elastic)  = a (intermediate  variable) 

D 

M „ 4 F „ 

“r  = + XB 

elastic  El  2 El 

2 2 
^ (I).  «,  F . I 

= F . A . zA  ^B 

— nr- 


Deflectlon  of  point  A due  to  vertical  F^  force 


Xa  = X3  + a • • cos 

‘h  ^A^'^zA  El  ^'*’a  ^A  • '’zA 


For  cos  ({1^  1.0 

*A  “ El  ‘•’a  6 ‘*’b^’^zA 


130 


SUMMARY  EQUATIONS 


For  un-bent  rod,  El  same  in  both  axes 


r xT 

11  1 0 0 

fF  ~1 

A 

3 6 

^A 

. 

_ 

1 i 0 0 

El 

rv  ^85 

1 

^A 

° ° 3 6 

F 

^A 

^.,51 

/b 

0 0 6 3 

F 

^B 

For  bent  rod,  with  point  A initially  at 

unstressed  position  having  x and  y components 

of  y^^  and  point  B having  corresponding 

initial  unstressed  deflection  components  of 

X-  , y write  out  deflection  equations 
D Ud 

in  algebraic  form. 


JG  ^^OA  “ ^ ^x. 


in-plane  bending 


torsion 


El 


2 ( OA  OB)  5 ^ 

2 £ 6 ■ 1 


Z-X  coupling 


Simi larly ; 


^A  ~ El 


^|f 

3 ^A  ^fi 


IG  ^^OA  - ^ ^OB^  ^v. 


!i3 

1 . ± 

( - 1 

^OA  ■ ^ObI  . 5 

^OB  1 

1 El 

2 ^ 

‘ ' 6 

a 

y 

\ 

1 
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For  Point  B (elastic  def  1 ctions) 


Note 


^ = ii  ( 5 F + 1 p ) 
B El  ^6  X.  3 X 
y A B 


f *0A  ^OB 

+ 1 

125 

El 

y 

1 ^ 

^ 3 

« 1) 

In-plane  bending 

(no  torsional  contribution 
at  B)  (1st  order  approx.) 


z-x  coupling 


^B  ~ El 


(■|  F + ^ F ) 
^ ^ ^B 


^nB 

El  4 Z 


+ - (— 1|F 
3 ' i / z. 


: Rod  bending  stiffnesses  EI^  and  EI^  are  designated  as  follows: 

EI^  is  stiffness  for  bending  rotations  about  x axis,  resulting  in 
y deflections 

Ely  is  stiffness  for  bending  rotations  about  y axis,  resulting  in 
X deflections. 
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FINAL  2-BEAM  GRAVITY  ROD  INFLUENCE  COEFFICIENT  MATRIX 


Numerical  Check 


For  manufacturing  tolerance,  min  radius  of  curvature  for  BI-STEM  > 1800.  foot 
per  NrS-2  spec. 


For  L = 2 i = 720" 

tie  deflection  = 12.0  in.  6 = 1.91021317° 

Add  24  in.  thermal  deflection  (per  R.G.  Langton  memo) 

Total  x^.  = 36.  in.  6 = 5.7248° 

OA 

R = 601.5  ft.  = 7218.  in.  (radius  of  curvature) 

XQg  = 9.005  inches  M = Yq^  = 0.0 

El  = 3000.  = El 
X y 

JG  =25. 

I = 360 . 

A 2-beam  Stardyne  finite  element  model  was  run  using  the  above  parameters, 
including  the  radius  of  curvature.  Two  straight  beam  element  connect  the 
base.  Node  B,  and  Node  A (top).  Run  TSNS  T6Z. 

An  additional  run,  using  twice  as  many  beams  and  nodes,  was  also  run  (Run 
TSNS  T9D)  . Both  results  are  presented  in  Tables  I-l,  1-2,  and  1-3. 


STARDYNE  CHECK  (,ASE,  GRAV.  ROD  INFLUENCE  MATRIX.  [E] 


TABLE  I- 10.  COMPUTED  [e]  ~ lbs /inch 


41,4  72. 

12,960. 

0 

0 

2073." 

rF  ^ 

^A 

1 

12,960. 

5 ,184. 

0 

0 

713. 

/ 

F 

"b  1 

^A 

> = 

0 

0 

46,138. 

12,960. 

0 

F 

^A 

0 

0 

12,960. 

5 ,184. 

0 

F 

^Xa^ 


TABLE  I- 11.  STARDYNE  2 BEAM  SOLUTION 
(INCLUDES  Z DEFLECTIONS) 


41,498. 
12  ,99  3. 
0 

0 

+2,465. 
+ 325. 


12,993. 


5 ,205.  0 


0 

0 

715. 

130. 


46,294. 

13,011. 

0 

0 


0 

0 

13,011. 

5,208. 

0 

0 


2464. 

715. 

0 

0 

149. 

18. 


F ^ 
^A 


V 


6x1 


5x1 
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TABLE  1-12.  Si'ARDYNE  4 BEAM  STATIC  CHECK  RESULTS 
FOR  FLEXIBLE  GRAVITY  GRAD.  ROD 

(COMPARE  TO  TABLE  2 & 1) 

(INaUDES  Z DEFLECTIONS) 


r 

41,494. 

12,992. 

0 

(?"■ 

0 

2,562? 

^A 

12,992. 

5,204. 

0 

0 

739. 

F 

^A 

0 

0 

48,336. 

13,886. 

0 

< 

F 

^A 

> 

^B 

0 

0 

13,886. 

5,354. 

0 

F 

^B 

^A 

■ 

2,562. 

739. 

0 

0 

162. 

F 

398. 

154. 

0 

0 

23. 

Run  TSidST9D 
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SHEARS  AND  MOMENTS  TRANSMIl > ED  TO  ROD  ATTACH  POINT  BY  FLEXIBLE  ROD  MODEL 


In  terms  of  previously  defined 
parameters  the  dynamic  loads  imposed 
on  spacecraft  by  rod  motions  are  as 
follows  : 

For  loads  in  Y-Z  plane: 

(shown  in  positive 
direction) 

Loads  on  spacecraft 
are  equal  and  opposite: 


V = F + F 

y ^A  ^B 

V = F 

z z 

M = - (?.  cos  (|),  + !!,  cos 
X A B 


i cos  • F 
B y. 


- <’'0A>  L. 


For  1st  order  non-linear  analysis,  update  moment  to  include  deflection  of 
point  A.  Assume  cos 


M = - 2 

X 


"A 

'a 


F - 8, 

Y/ 


'y^  - ^^OA  - 


L 


I 

I 

f 
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Loads  in  X-Z  Plane 


I 


Loads  shown  at  base  of  rod  to 
maintain  equilibrium. 

Loads  on  spacecraft  are 
opposite  and  equal. 


V = F + F 

* ’'a  ’'b 

(previously  calculated  = F^ 

= +(«.  cos  41.  + cos  <)>,,) 
y A c X 


+ i cos  <j>T.  F 


* L<-0A> 


For  1st  order  non-linear  analysis,  update  last  term  to  include  current 
deflection  of  point  A (tip  mass)  which  could  be  changing  fastest.  Also 
1 


cos  (b,  = cos  (ti„ 
A B 


M=2JIF  +ilF  +F  (x-.  - xj 
y x^  Xg  2 OA  A 


Moments  about  Z Axis 


From  top , looking  down 


M = y~.  . F + y«„  • F 
z -^OA  X,  Ob  X 
A o 


OA  y. 


^B  ^ 


XA 


For  non-linear  effects,  qpdate  total  distance 
from  Z axis . 


X 


oA 


! 

I 

I 

I 

I 

\ 

I 
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Non-linear 


^^OA  - V ^.  ■ 


x„)F 

y. 


In  matrix  notation,  loads  on  spacecraft  at  rod  attach  point  are: 


r 'N 

V 

1 

1 

0 

0 

0 

0 

F 

X 

V 

0 

0 

1 

1 

0 

0 

^A 

F 

y 

V 

0 

0 

0 

0 

1 

1 

^B 

F 

z 

M 

X 

>■ 

0 

0 

-n 

-!!, 

(yoA-^A^ 

(yoB-yfi^ 

yA 

F 

yB 

M 

y 

2^ 

i 

0 

0 

F 

^A 

M 

z 

y-*S 

o 

> 

> 

^yoB-^B^ 

-^^0A-^A> 

0 

0 

F 

V.  ^ 

L_ 

— 

where 


X, , x_  , y,  , y„  are  updated  rod  elastic  deflection. 

A D A D 

For  non-linear  Lagrangian  analysis  or  do  not  include  for 
linear  (small  deflection)  analysis. 
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NUMERICAL  CHECK,  REACTIONS  AND  MOMENTS 


Linear  analysis 

CALCULATED  VALUES,  STARDYNE  CHECK  CASE 


r 

V 

X 

1.0 

1.0 

0 

0 

0 

0 

^A 

V 

y 

0 

0 

1.0 

1.0 

0 

0 

F 

V 

z 

> = 

0 

0 

0 

0 

1 

1 

F 

^A 

M 

X 

0 

0 

-720 

-360 

0 

0 

F 

^B 

M 

y 

720 

360 

0 

0 

36 

9 

F 

"a 

M 

0 

0 

-36 

-9 

0 

0 

F 

— 

— 

COMPUTED 

reactions 

FROM  STARDYNE 

2 BEAM 

model  (run  USNS  T6Z) 

r '\ 

V 

1.00 

1.00 

0.0 

0.0 

0.0 

0.0 

A 

^A 

V 

y 

0.0 

0.0 

1.00 

1.00 

0.0 

0.0 

F 

V 

z 

Y- 

0.0 

0.0 

0.0 

0.0 

1.0 

1.0 

1 

F 

^A 

M 

X 

0.0 

0.0 

-720.00 

-360.45 

0.0 

0.0 

\ 

F 

^B 

M 

y 

720.0 

360 . 45 

0.0 

0.0 

36.0 

9.006 

F 

"a 

M 

0.0 

0.0 

-36.00 

-9.006 

0.0 

0.0 

F 

L^bJ 

lAO 

1 


1-3  DISCRETE  DAMPING  MATRIX  DERIVATION  AND  CAI  CULATION 


I 

( 

I 


Convert  Modal  Damping  to  Discrete  Damping 


Matrix  equations  of  motion  for  free  vibration: 


[m]  {x}  + [c]  fx)  + [k]  {x)  = {0}  (I-l) 

Let  {x}  = [(()]  {q}  (I-?.) 

where 

[(fl  = set  of  mode  shapes 
{q}  = modal  displacements 

Also  {x}  = [())]  {ql  (1-3) 

(x)  = [(})]  {q-l  (I-.-i 


Substituting  Eqs.  1-2,  1-3,  and  1-4  into  Eq . I-J  and  premultiplying  by  [<() 
the  result  is 

[Mj,q]  {q}  + [Mg^  2 Sto^llq}  + [M^,^  fq}  = 10}  (I-S) 


where 

[%]  = [M]  [^] 


(Generalized  Mass  Matrix)  or  equivalent 
(Ref.  STARDYNE  Theoretical  Manual, 
DYNREl  Analytics) 


For  solution  by  the  normal  mode  method,  Eq . 1-5,  the  damping  must  be  assume^, 
to  be  proportional  to  the  mass  or  to  the  stiffness.  For  the  first, 

[C]  = 2 Dl]  [m]  (l_r,) 


If  so,  then 

uf  [c]£^]  = 2 [B] 


(1-7) 
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Note : 

(Corollary  for  one  d.o.f.,  oscillator) 

Ref:  Thomson,  W.T.,  "Vibration  Theory  and  Applications,"  Prentice-Hall, 
1965,  p.  51-55. 

= 2 m 10^  (critical  damping  coeff.) 

C ’ 

P = — (damping  factor) , 

c 

C'=2pmw  (8) 

n 

From  Equation  1-7,  assume  a value  of  ,005  for  P,,3„,  .,,,  3 . 

T — 1 12  n 

Premultiply  both  sides  of  1-7  by  . (Inverse  at  transpose  of  [(j)]). 

[I]  [C]  [<^]  = [lif  2 [S]  [m^q]  (9) 

Post-multlply  both  sides  of  Eq,  1-9  by 

[C]  = [2  8][MgQ]  [4.]"^  (10) 


Evaluate  numerically. 

I 


Mode  No. , n 

lbs. 


H 

z 


1 

2 


M = 


Uf  = 


[2  S]  = 


7.7969  1.000 

.47307  -.01526 

T.OOOO  -.01526 

_.  31352  1.0000_ 

1.0000  .31352 

-.01526  1.0000_ 

~.01  0 
_ ^ 

“7.  7969  0 

0 .47307 


.31352  .00550005 

1.000  .13555 

Columns  are  for  a 

mode  (11) 

(12) 


lbs . 


[2  6]  [M^q]= 


.01 

0 

7.7969 

0 

0 

.01 

0 

.47307 

.077969 

0 

0 

.0047307 

1 

1.000 

.31352 

-1 

1 

1.000 

- 31352 

-.01526 

1.000 

+.01526 

1.000 

(13) 
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whe  re 


de  termin  u;  t 


^11^22  ~ ^12  ■ ^21 
= (1.0)  + (.31352) (.01526) 

= 1.004  7843 


.995238 


.0151873 


-.312027 


.995238 


1.000 


.01526 


.31352  1.000 


1.000 


.01526 


-.31352  1.000 


|A|  = (1.0)^  + (.01526)(. 31352) 


1.0047843 


995238 


-.312027 


.0151873 


.995238 


Performing  numerical  multiplication,  from  Eqs.  1-13  and  1-14 


[Uf]  [2B][M£p] 


.995238 


.0151873 


0775977 


.0011841 


-.312027 


.995238 


-.001476 


.00470817 


.077969 


.0047307 


Then  post-multiplying  by  Eq . 1-15, 


.0775977  -.001476 


0011841 


.995238  .0151873 


.00470817  -.312027  .995238 
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.07768884  -.00029057 
-.00029057  .00470374 


386.4 
(Divide  by  386.4) 


ib -sec/inch 


For  1 d.o.f.  system,  f = .0055  H. 

n z 


u = 2'  f = .0346  rad/sec. 
n n 


Using  Eq . 1-8 

C = 2 e = 2 (.005)  (7.725/386.4)  (.0346) 

= 6.917  X 10  ^ Ib-sec^/in.  — 

sec 


If  divided  between  Node  A and  Node  B 

C.  = C„  = 3.4586  X 10"^ 

A B 


Dividing  coupled  daiiping  matrix  by  386.4  (since  Stardyne  is  expressed  in 

EQ 


lbs.) 


[c]  = 


2.01058  X 10 
-.75199  X 10' 


-.75199 

X 10”^ 

1.21732 

X 10”^ 

Ib-sec/in . 


Note:  if  single  d.o.f.  damping  were  based  on  second  mode  only 
C = 2 B mg  ^2  (^2  = 2tt  (.  13555) 


C = (.01) 


.45 

386.4 


271  (.  13555) 


= 9.9187  X 10 


-6 


(see  below) 
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For  single  d.o.f.  system 


Damping  factor  = C c/in 

Damping  force  = C 
For  sine  vibration 
* = uX  in/sec 

Damping  force  = C w X Ib-sec/in  • rad/sec  in.  = lbs, 

C'  = 2 6 m = 2 X (.005)  = 1,9992  x lO"'^ 

1,1  A 386.4 


For  large  mass . 


C ’ 

^2,2 


2 6 


= 2 X .005  X 


.45 

386.4 


1.1646  X 10 


-5 


CHECK  damping  MATRIX 


Assume  constraiined  vibration  mode,  at  resonance, 
tip  anplitude  = 1.0  sin  w t 


X,  = 1.0 
A 


X = w X.  = .03455 
A A 


.31352 


= u x„  = .010834 

..2  ® 2 ® 

X.  = -0)  X.  = -.0011942  in/sec 
A . A _ 

Xg  = -0)^  Xg  = -.00037441  in/sec 


f = .0055  H 

z 

0)  = 271  f = .03455 


Inertia  forces; 


M,  0 

A 


0 M., 


7.725 

386.4 

.45 

386.4 


X ,0011942 


X .00037441 


-2.3875  X 10 


-4.3604  X 10 


lbs, 


I 
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For  stiffness  forces,  must  invert  2x2  Influence  matrix 


S tl f fness  forces  : 
Evaluate  numerically: 


1.10225  X lO"^  -2.7557  x 10“^ 

[k]  = lbs /inch 

-2.7557  x 10"  8.8183  x 10"_ 
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NODE 

INERTIA,  lbs. 

STIFFNESS.  Ib. 

DAMPING 

-5 

-5 

^-6 

A 

-2.3875  X 10 

-7 

2.38315  X 10 

-7 

6.9384  X 10  ° 
-7 

B 

-4.3604  X 10 

8.9947  X 10  ' 

1.059  X 10 

At  resonance. 

inertia  + stiffness 

forces  are  equal  and 

opposite . 

1 

1 

I 
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APPENDIX  J 


DETAILS  OF  '.OLAR  ARRAY  FLEXIBILITY  MODELING 

SOLAR  ARRAY  TORSION 

The  solar  array  torsional  equations  of  motion  are  described  by  a solution  of 
the  following  system  of  matrix  equations 

[m]  {e}  + [C]  {6}  + [k]  {0}  «=  {T(t)}  (J-1) 

where  the  vector  {0}  is  the  8^  rotation  of  each  solar  panel  and  the  drive 
motor. 

Applied  torques  T(t)  may  be  zero  at  all  locations  except  the  drive  motor. 

The  above  equations  in  discrete  coordinates  are  coupled  becaxise  the  stiffness 
matrix  [k]  is  coupled.  The  damping  matrix  [c]  may  also  be  coupled.  For 
six  solar  panels  plus  the  drive  motor,  the  equation  has  at  least  seven 
degrees  of  freedom.  Now  malce  the  following  transformation: 

{0}  = [()>]  {q}  (J-2) 

where  the  columns  of  [<{>3  are  the  eigen  vectors  which  result  from  the  STARDYNE 
solution  of  the  relationship 

[K]  {<^}  = J [m]  {<!>} 

Substituting  equation  J-2  and  its  derivatives  into  Eqn . j_i  and  premultiplying 
by  the  transpose  of  [(j)]  results  in  the  followiuj'.  set  of  uncoupled  equations  in 
terms  of  the  generalized  coordinates  {q}: 
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•T-i;) 


[Meq]  (q}  + (q } + = [41]^  (T(t)} 

Note  that  by  orthogonality  relationships,  the  oif-diagonal  terms  of  the 
teneralized  mass,  damping  and  stiffness  matrices  are  zero. 

The  STARDYNE  solution  for  the  NTS-2  model  with  flexible  solar  arrays  and 
flexible  gravity  gradient  rods  has  six  rigid  body  modes  and  one  mechanism  type 
rigid  body  mode  because  the  solar  arrays  are  considered  free-wheeling  in  6^. 
The  "free-wheeling"  solar  array  mode  in  the  STARDYNE  modes  is  included  as  a 
linear  combination  of  all  seven  STARDYNE  rigid-body  modes. 


For  simulation  purposes , the  free-wheeling  solar  array  mode  will  be  replaced 
by  a single  equivalent  rigid  body  6^  mode  whose  anplitude  is  normalized  to 
1.0  and  whose  generalized  mass  is  thus  equal  to  its  total  inertia  about  the 
solar  array  9^  axis. 


The  elastic  modes  are  taken  from  the  STARDYNE  run  TSNSTRL.  The  only 
significant  solar  array  elastic  modes  are  the  symmetric  torsional  modes  at 


1.095  H and  1.  89  8 H . 
z z 

Thus  the  total  number  of  generalize 
defined  as  follows: 


Generalized  Mass 

[M^q]  = Uf  [M]  M 


coordinates  to  be  solved  is  three, 

= rigid  body  torsion  mode 
= first  symm.  torsion  mode 
q^  = second  symm.  torsion  mode 
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The  diagonals  of  the  generalized  weight  matrix  are  printed  by  STARDYNE  and 
are  claculated  for  modes  normalized  to  a maximum  displacement  or  rotation  of 
unity.  For  the  above  three  solution  modes: 


MODE 

NO. 

equiv.  weight 

M , Lbs . 

eq 

EQUIV.  MASS 

M , Ib-sec^/ft 
er 

FREQUENCY 

H 

z 

0) 

n 

rad/sec 

w 2 

^ 2 
(rad/sec) 

1 

13956.7 

3.01 

0 

0 

0 

2 

9284. 

2.002 

1.095 

6.880 

47.  339 

3 

69  71. 

1.5034 

1.898 

11.925 

142.2 

Premultiplying  Eqn . J-3  by  [W^q]  the  response  equation  becomes 

T 

[I]  {q}  + [2pa)J  {q}  + [0)^2]  {q}  = ^ 


Extracting  only  the  0^  rotations  at  the  drive  motor  (node  20)  , the  -Y  outer 
panel  center  (node  2) , and  the  +Y  outer  panel  center  (node  31) , these 
rotations  are  related  to  the  generalized  responses  by 


1 


9 drive  motor 
6 -Y  tip  mass  ) 

y 

6 +Y  tip  mass 

y 


Equation  J-4  in  numerical  form  reduces  to  (p  = 


O-L) 


1.0 


1.0 


1.0 


r‘ 


0.0 


.0688 


.11925 


0.0 


47.334 


142.2 


.33222 

-.4995 

-.3326 


(T,  . (t)) 

drive 

motor 


(•j-r.) 
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Equation  J-6  can  be  solved  as  three  uncoupled  cond  order  linear  dif- 
ferential equations. 


TABLE  J-1.  -Y  SOLAR  ARRAY,  (BENDING  ONLY) 


STIFFNESS  MATRIX 

(Units  = lbs. /ft.)  For  motion  relative  to  drive  motor 


[k]  {x} 


.0085617 

-.004860 


.004860' 

.630218 

where 

= motion  of  center  of  outer  panel,  normal  to  plane  of  panel 
= motion  of  center  of  outer  panel,  in  plane  of  the  panel 


MASS  MATRIX 

2 

(Units  = slugs  (lb-sec  /ft) 


[m]  {R}  = 


.42176 

0 

0 

.52608 

[hj 

DAMPING  MATRIX 


[C]  (X)  = (^) 


P = .005  (%  of  1 percent  crticial  damping) 
K 


.0006009  0 

0 .005758 


(Units  of  [c]  are  Ib-sec/ft.) 


TABLE  J-2.  +Y  SOLAR  ARRAY,  (BiLDING  ONLY) 


STIFFNESS  MATRIX 
(Units  = lbs. /ft.) 

[k]  {x}  = 

MASS  MATRIX 

2 

(Units  = slugs  (lb-sec  /ft)) 

,42176  0 


[m]  f x)  = 


For  motion  relative  to  drive  motor 


.00  86  70 

-.010  705 

r»2) 

-.010705 

1.37791 

Ui 

0 


.52608 


DAMPING  MATRIX 
(Units  = (Ib-sec/ft) 

[c]  {xl  = 


.0006047 

0 

0 

.008514 

1^2  J 

TABLE  J-3.  NTS-2  SOLAR  ARRA  ■< 


COMPARISON  OF  MAXIMUM  OUTER  PANEL  DEFLECTIONS 
FOR  STEP  INPUT  YAW  TORQUE  OF  .04  FT-LBS  FOR 
3 SECONDS,  RESPONSES  ARE  NORM.  TO  PANEL,  AT  CENTER 


DYNREI 

RUN  NO. 

MAX.  REL. 
RESPONSE,  NODE  2 
INCHES 

TIME  OF 

MAX.  RESPONSE, 
SECONDS 

MODE  NOS. 
INCLUDED 

IN  SOLUTION 

COMMENTS 

TSNSTCO 

-.009843 

1.40 

10,  14,  17,  24 

(Tape  X3350) 

Includes  1st 

and  2nd  anti- 
symmetric 

array  bending 

modes 

TSNSTIl 

-.009855 

1.40 

10,  14,  17 

(Tape  X3350) 

2nd  anti- 
symmetric 

array  mode 

de leted 

TSNSTJO 

-.009916 

1.40 

10,  14,  17 

Simplified 

bending 

model 

i 

1 
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Figure  J-1.  Stardyne  Math  Model  for  Flexible  Solar  Arrays 


Figure  J-2. 


Gravit' 


Gradient  Rod  and 


56 


'.|-acecraft  Definition 


1 


Figure  J-3.  ^olar  Array  Bending  Flexure  Model 
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REL  DISP  REL  VEL  REE  AfX  nz 


0 2.0000  4.0000  6.0000  8.0000  10.0000 


0 2.0000  ^.0000  6.0000  8.0000  10.0000 

TIME  (SECOtIDS) 


Figure  J-4.  Dynamic  Resi'onse  of  Outer  Panel  (Node  2)  to  Step  Yaw 
Torque  Input  of  .04  ft-lbs  for  3 seconds 


158 


T 

I 

I 

1 


15.10 


Figure  J-5.  Solar  Array  Drive  Motor 
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REL  DISP  REL  VEL  REL  ACC  GZ 


I 

i 


6.0000x10 

4.0000x10 
2.0000x10 
0 

-2.0000x10 

-4.0000x10 

-6.0000x10 


-01 


-01 

-01 


-01 

■01 

-01 


6.0000x10 

4.0000x10 

2.0000x10 

0 

-2.0000x10' 

-4.0000x10 

-6.0000x10 


-02 

-02 

-02 


1.00000  2.00000  3.00000  4.00000  5.00000 


■02 

-02 

-02 


1.5000x10 
1.2000x10 
9.0000x10 
6.0000x10 
3.0000x10 
0 


Figure  J-8.  Solar  Array  Torsion  Model , Mode  2 = Outer  Panel 
Dlsp.,  Vel.,  & Accel.  Response  in  Radian  Units 
to  Four-Pulse  Excitation 
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REL  DISP  REL  VEL  REL  ACC  GZ 


REL  DISP  REL  VEL  REL  ACC  GZ 


IIRIVE  MOTOR 
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igure  J-12.  Mode  Plots — Torsion  (Pirst  Order  Modes) 


DRIVE  MOTOR  (ZERO  I 


7 'dsia  X nvaow 


(H31JI/S3HdlJl) 

•dsiu  X nvaow 


TABLE  J-5.  STARDYNE  MODAL  DATA:  FLEXIBLE  RODS— 2 BEAMS/ROD 

FLEXIBLE  ARRAYS~4  DOF/ARRAY  (EACH  ARRAY) 


rray 


To  kill  of  .60  rad/sec  in  325  msec 


6 = 0 - e.t 

o 

T 

0 = .60  - (0.325) 

Tx.  =•  “ 1.846154  In-lb. 

b • jZj 


^ -1.846154 

T “ 12 


-.1538462 


NTS-2  SOLAR  ARR.iY  STRUCTURAL  ANALYSIS 


o A linear  Open  Loop  mode]  of  total  NTS-2  spacecraft  was  run 
on  STARDYNE, 

o All  sub-structure  stiffnesses  combined. 

o Contiined  with  previous  defined  beam  models 
(2  beams/rod)  of  gravity  gradient  rods. 

o Above  model  was  reduced  to  varying  degrees  of  complexity  by 
selecting  fewer  lumped  inass  points  f solar  arrays. 

RESULTS 

o Solar  array  flexible  torsion  model  should  have 
4 D.O.E.  + drive  motor  V. 

o Solar  array  flex,  bending  model  can  have 
2 DDK  per  array  with  good  results, 

4 DDK  per  array  with  better  results. 

o Open- Loop  transient  response  rins  to  be  made 
to  guide  decision  i*i  bending  model. 
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NTS-2  SOLAR  ARRAY  STRUCTURAL  ANALYSIS 


FOR  EACH  ARRAY,  CONSIDERED 

o 3 Solar  Panels  (identical) 

o 1 Interface  Panel 

o 4 Hinges  (/^sumed  Identical) 

o Drive  Tube  (L  = 35.53  in.,  -Y  Array) 

(L  = 14.87  in.,  +Y  Array) 

Each  identical  structure  was  modeled  in  detail  by  finite  element 
analysis  (STARDYNE) , and  reduced  to  a sub-structure  stiffness 
matrix  element. 

Sii)structures  for  -Y  array  checked  out  for  comparison  with  NRI. 
vibration  test 

o When  gravity  field  stiffness  included,  natural 
frequency  in  X and  Y agreed  well. 

o Model  is  somewhat  stiffer  in  bending  in  Z direction 
(plane  of  panel)  due  to  lack  of  defn.  of  brkt.  at  end 
of  drive  tii>e. 


/ 

i 
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PITCH  SENSOR  YAW  SENSOR 


Figure  J-17.  Deployed  Solar  Array  Configuration  (NTS-2 
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X2  VIEW  ROTATED  AXES 


5.05 


STARDYNE  MODEL 
NTS-2  HINGE 


Figure  J-19.  Stardyne  Model  NTS 
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TABLE  J-6.  SUMMARY  OF  SOLAR  ARRAY  CANTILEVER  VIBRATION 
TEST  RESULTS  @ NRL  VS.  STARDYNE  MODEL 
(GRAVITY  VECTOR  ALONG  Y AXIS) 


1ST  TORSION 
FREQ. 

0 

7 

H 

z 

p 

1ST  X BENDING 
FREQ . , NORMAL 
TO  PANELS 

H 

z 

1ST  Z 

BENDING 

FREQ. 

H 

z 

1 

INFORMAL  NRL 

TEST  DATA,  (REF. 

LARRY  TURNER) 

.40 

.44 

00 

STARDYNE , 

ORIG.  HINGE  STIFF. 

I G FIELD 

.408 

.48 

2.07 

STARDYNE , 

MODIFIED  HINGE  STIFF. 

I G FIELD 

. 384 

.47 

2.05 

STARDYNE, 

ORIG.  HINGE  STIFF 
"ZERO  G" 

.408 

.28 

2.04 

STARDYNE, 

MODIFIED  HINGE 
"ZERO  G" 

1 1 

1 

.384 

.26 
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TABLE  J-8.  STAPDYNE  MODELS,  TOTAL  NTS-2  SPACECRAFT 
EFFECTS  OF  COMPLEXITY  OF  SOLAR  ARRAY 
MODEL  SECOND  ORDER  MODES  COMPARISON 
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TABLE  J-9.  STARDYNE  MODAL  DATA  (concluded) 


TABLE  J-10.  STARDYNE  MODAL  DATA:  FLEXIBLE  RODS — 2 BEAMS/ROD 

FLEXIBLE  ARRAYS~3  DOF  AT  INNER,  OUTER  PANELS 


-112.23  in. 


6Y3  rotation,  inner  panel,  +Y 

0Y4  ROTATION,  OUTER  PANEL,  +Y 


Cx[K1{0}=CxKin  FT-LB/RADIAN) 

where  C = 1/12 


211.065 

" evi" 

(SYMMETRIC) 

-211.131 

352.756 

1 

eY2 

1 

12 

0 

-141.625 

283.281 

6YM 

0 

0 

-141.788 

352.919 

eY3 

0 

0 

0 

-211.131 

211.065 

1 

. eY4  . 

1 

(Final  model  is  a modal  solution  using  the  more  detailed 
structural  model  for  torsion.) 


Figure  J-21.  Solar  Array  Torsion  Flexure  Model 
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APPENDIX  K 


MDDIFIED  NTS-2  SOLAR  ARRAY  TO', ION  MODEL 

Hinge  stiffness  increased  from  326.  to  92,000.  in- lb /rad.  Drive  tube  O.D. 
Increased  from  1.25  to  1.75  Inches.  Drive  tube  wall  thickness  increased 
from  .049  to  .063  inches.  (Above  information  from  Larry  Turner,  NRL  by 
telephone  calls  1/18/77  and  1/26/77.) 

Previous  model  of  solar  arrays  solved  for  torsion  only  (degree  of  freedom  5 
= GY). 


For  existing  S&R  simulation  model,  update  as  follows: 


Mode  No. 
(Symmetric  Only) 

Equlv.  Weight 

M Ib-in'^ 

eq 

M 

eq 

Equiv.  Mass 
Slug-ft. 

Freq  . 

H 

z 

CO 

n 

rad/sec 

2 

Ul 

2"  2 
rad  /sec 

139  4 4. 

3.007 

0 

0 

0 

7152.5 

1.5426 

7.9635 

50.036 

2503.6 

6289.1 

1.3563 

15.5315 

97.587 

9523.3 

The  differential  equation  of  motion  of  the  solar  array  drive  system  (open  loo 

loop)  is  given  in  terms  of  the  generalized  coordinates,  q,  as: 

T 

I (q)  + 2 po)  fq)  + w ^ fq}  = 

n ^ n ^ M 

eq 

For  0 = 0.005  (h  percent  critical  damping) , the  generalized  equation  reduces 
to 
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1.0 

3.007 

i 

-1.00 

1.5425^ 

-.8856 

1.3563 

^ . 33255 


UT 


drive  motor 


(t))  = < -.64830 

..-.65295  V 


' ^^drive  motor 


This  replaces  Eq.  6 in  previous  memo.  The  discrete  responses  of  several 
locations  of  interest  are  taken  from  the  mode  shapes  and  given  below: 


LOCATI(»I 

STARDYNE 

6 Rotation 

y 

(Radians)  per 

Unit  q 

MODE  NO. 

Mode  1 
(Rigid  Body) 

Mode  2 

Mode  3 

Drive  Motor 

20 

1.0 

-1.0 

-.8856 

Tip  of  -Y  outer 
panel 

1 

1.0 

1.0 

-.8856 

Pitch  sensor 

3 

1.0 

.5683 

.32105 

Tip  of  +Y  outer 
panel 

32 

1.0 

1.0 

-.8856 
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TABLE  K-2.  ANTl-SV^IMETRIC  TORSION  MODE 
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TABLE  K-3.  SY>D-1ETR1C  TORCIOH  MODE 
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TABLE  K-A.  2Nn  ANTI- SYMMETRIC  TORSION  MODE 
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TABLE  K-5.  2ND  SYMMETRIC  TORSION  MODE 
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APPENDIX  L 


SOLAR  ARRAY  TORSION  MODES  , 

SOLUTION  FOR  RESPONSE  TO  TORQUE 
MOTOR  ANGULAR  ACCELERATIONS 

Summary  input  data 

Based  on  updated  solar  array  hinge  torsional  stiffness,  the  arrays  are 
assumed  to  have  symmetric  responses  for  a rigid  torque  motor  connection  to 
the  spacecraft. 

The  modes  of  interest  are  obtained  from  STARDYNE  finite  element  structural 
model,  using  only  the  torsional  masses  and  spring  stiffness  data,  including 
updated  hinge  stiffness.  For  a locked-up  or  grounded  torque  motor,  the 
following  data  is  obtained  for  the  two  modes  of  interest: 


Mode  No . 
(STARDYNE) 

Mode  No. 

(S  & R) 

Equiv.  Wei^t 

M lb- in 

eq 

M 

eq 

Equiv.  Mass 
Slug-ft. 

Freq  . 

H 

z 

0) 

n 

rad/sec 

— 

OJ  2 

2^  2 
rad  /sec 

2 

1 

100  72.29  7 

2.17225 

2.5908 

16.2784 

264.987 

4 

2 

7719.405 

1.6648 

9.0358 

56.774 

3223.281 

The  differential  equations  of  motion  of  the  motor  array  drive  system  in 
terms  of  the  generalized  coordinates  q is  given  by: 


[ i]{q}  + [2 

where 


(L-1) 


e 

b" 


angular  acceleration  of  drive  motor  (absolute) 
I Torques  applied  to  drive  motor,  not  ^ ^ 

pitcL’ 


spacecraft /pitch 


(L-2) 

acceleration 
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Note  that  if  torque  motor  is  moving  and  torque  is  applied  to  cause  + S' 

3 r r sy  s 

then  a negative  torque  is  applied  to  the  spacecraft. 


6 

sc 


pitch  moment  of  inertia  of  spacecraft,  without  arrays 

-(f  Torques  applied  to  drive  motor)  + all  other  pitch  torques 

I 


Also  note  that  the  net  torque  applied  to  drive  motor  = T u j ~ Tt,  i ^ „ 

^ applied  brake  torque 


■■  Drive  motor  Torque  , " 

3 = 1 + 

array 


(L-3) 


+ 


.16278  0 

0 .56774 


264.987  0 

0 3223.28 


-.  80  76 
-.2076 


(+  e^) 


pitch  sensor 
T 

e = -S  + e 

A I sc 


= .95183  q^  - .46754  q^ 


For  P = 0.005  (h  percent  critical  damping)  the  generalized  torsional  equations 
of  motion  become 
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The  discrete  responses  of  several  locations  of  Interest  on  the  solar  panels 
are  taken  from  the  mode  shapes  and  the  generalized  coordinate  responses 

{6}  = {q}  -H  6 ^ (L-5) 

where  B = absolute  angular  dlspl.  of  drive  motor  / 


STARDYNE 
MODE  NO. 

LOCATION 

<j)  MATRIX  TERMS 

20 

Drive  motor 

0.000 

0.000 

1 

Tip  of  -Y  outer  panel 

1.000 

-.98418 

3 

Pitch  sensor 

.95183 

-.46 754 

32 

Tip  of  +Y  outer  panel 

1.000 

-.98418 

for  both  sides 


APPENDIX  M 


COMBINED  EARTH  SENSOR  MODEL 


INTRODUCTION 

This  appendix  describes  the  model  developed  for  the  GPS  combined  earth  sensor. 
The  model  can  be  used  to  evaluate  acquisition  as  well  as  normal  operation. 

This  model  is  based  on  information  obtained  from  Barnes  Engineering  Company 
at  a meeting  in  St.  Petersburg,  Florida  on  March  13,  1975,  and  a subsequent 
phone  call,  and  the  Rockwell  procurement  specification  for  the  sensor. 

SENSOR  DESCRIPTION 

The  Combined  Earth  Sensor  is  a two- axis , static  sensor  which  incorporates 
two  horizon  crossing  indicators  for  spinning  mode  attitude  sensing.  The 
adjective  static  means  the  sensor  uses  no  moving  parts  in  the  measurement 
of  attitude. 

The  two-axis  sensor  consists  of  an  optical  system,  24  thermopile  detectors, 
and  processing  electronics.  The  24  detectors  are  arranged  on  the  focal 
surface  of  the  optical  system.  The  earth  image  is  tangent  to  the  inner 
edges  of  the  A detectors  at  null  at  the  design  altitude.  Two  roll  and 
two  pitch  measurements  are  computed  sequentially  by  switching  sums  and 
differences  of  various  detector  outputs  to  the  ininit  of  an  integrator. 

Each  measurement  is  sampled  and  held  for  comparison  check  and  possible 
output.  A new  roll  and  new  pitch  measurement  is  computed  each  1/4  second. 

Roll  and  pitch  computations  are  separated  by  1/'  second.  The  S detectors 
I are  for  space  reference  and  are  used  with  the  etectors  for  accurate 

I 
I 
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measurements  over  +1^  In  the  track  mode.  The  nd  B detectors  are  used 
in  the  acquisition  mode  for  +4*^  range.  The  th  ir.opile  detectors  have  a 
time  constant  of  1.4  seconds  and  are  sensitive  lo  the  14  to  16  micron 
wavelengths  (CO2  band).  The  sensor  continuously  tests  all  signals  and 
rejects  spurious  or  incorrect  signals  due  to  sun,  moon,  or  other  anomalies. 

The  two  horizon  crossing  Indicators  consist  of  3 degree  by  3 degree  square 
fields  of  view.  One  diagonal  of  each  square  is  parallel  to  the  scan  path. 
One  field  is  10°  from  the  -X  axis  toward  +Z  and  the  other  is  10°  from  -X 
toward  -Z  (Spin  rotation  is  about  the  Z axis).  The  detectors  are  pyro- 
electric with  insignificant  time  constant. 

MODEL  DESCRIPTION 


The  model  consists  of  an  initialization  and  two  normal  operation  sections. 

The  initialization  section  performs  one-time  computations  and  clears  outputs 
and  counters.  One  normal  operation  section  is  for  the  horizon  crossing 
indicators  used  in  the  spin  phase.  The  other  normal  operation  section  is 
for  two-axis  earth  sensing. 

Characteristics  of  the  horizon  crossing  indicator  model  are  as  follows: 

• Inputs  to  the  model  are  spacecraft  position  vector  from  center 

of  earth,  and  spacecraft  attitude  direction  cosine  matrix.  Both 
inputs  must  be  relative  to  the  same  reference  frame. 

• When  a sensor  line  of  sight  crosses  ttie  edge  of  the  earth,  the 
time  from  the  crossing  to  the  sample  time  is  output.  I f no 
crossing  occurs  in  the  sample  interval,  ttie  output  is  set  to 
-At. 

• If  a sensor  line  of  sight  crosses  the  earth  very  close  to  a tangent, 
there  may  not  be  an  output.  The  sample,  '.t,  should  be  small 
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enough  so  that  this  loss  of  output  does  not  happen  when  attitude 
determination  is  critical.  Select  At  as  follows: 

^2 


At 


< JL^/i  . ^ 
- W V r 


tan 


-1 


[(1  - 

e 


1] 


where 

W 


RPM 


= spin  rate  in  rad/sec 

= ^ W 
60  RPM 

= spin  rate  in  rev/min 

= satellite  altitude  in  n.  mi . f 100  h £ 10897) 

= earth  radius  in  n.  mi.  (3440) 

= distance  from  earth  center  to  chord  of  detector 
Path  in  n.  mi. 


e.g 

let 

W 


RPM 


= 100 
= 10897 
= 3440 

= 3405.6  = 0.99  r 


then  At  < 0.00326  seconds. 


o The  equations  used  in  the  flow  diagram,  are  derived  below 


' R X LOS  1 

SIN  (ANG2) 

R • LOS 

~~wr  = 

-COS  (ANG2) 

ANG2  = TAN~^ 

SIN  (ANG2) 

COS  (ANG2) 

AKT/^O  — . T'AM  ^ 

|R  X LOS| 

ANUz  = lAW 

-IT  • LOS 

20A 


ANGl  = SIN 


= TAN 


-1 


r 

e 


r^/lRi 


[1  - 


= TAN 


[|^ 

ANG3  = ANG2  - ANGl 


2 _ , 23'. 

e 


Note  that  3440  + 100  ^ |R,  10897  + 3440  n.  mi. 

13.8°  ^ ANGl  £ 76.4°. 

Characteristics  of  the  two- axis  earth  sensor  model  are  as  follows: 

• The  sensor  output  is  updated  every  1/8  second. 

• The  sensor  model  has  four  states  whicli  are: 

State  1 Compute  Roll  1,  Output  Pitch  2 if  valid 

State  2 Compute  Pitch  1,  Output  Roll  1 if  valid 

State  3 Compute  Roll  2,  Output  Pitch  1 if  valid 

State  4 Compute  Pitch  2,  Output  Roll  2 if  valid 

The  state  is  advanced  every  1/8  second. 

Each  detector  area  is  divided  into  1°  x 1°  squares.  The  area 
coverage  of  each  1°  x 1°  square  by  the  earth  is  computed  from  the 
distance  between  center  of  the  earth  circle  on  the  focal  surface 
and  the  center  of  the  square.  The  orientation  of  the  square  is 
not  considered. 

The  position  of  the  earth  image  on  the  focal  surface  (B&C)  is 
computed  from  true  attitude.  The  diameter  of  tiie  earth  image 
on  the  focal  surface  (RE)  is  computed  from  satellite  altitude. 
Measurements  on  the  focal  surface  are  in  equivalent  radians. 
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• Each  detector  output  Is  the  sum  of  tli  coverages  of  1°  x 1° 
squares  within  the  detector. 

• The  sensor  output  is  computed  from  the  detector  outputs  according 
to  the  functions  in  the  electronics. 


• Functions  of  the  system  programner  to  reject  spurious  or  Incorrect 
signals  is  Included. 

• Subroutine  must  be  Initialized  by  calling  with  MODE  = 0 prior 
to  operating  in  two- axis  mode. 

• The  time  increment.  At,  should  be  selected  from  one  of  the 
following  to  synchronize  with  the  sensor  update  frequency: 


Recommended  At 
0. 125 
0.0625 
0.025 
0.0125 
0.01 
<0.01 


Increments  per  Update 
1 
2 
5 

10 

12.5  I Synchronization  error 

>12  5 ( sec.  is 

— ‘ J insignificant 


EARTH  RADIANCE  VARIATION 

The  variation  of  radiance  over  earth  surface  for  the  14-16  CO^  band  has 
been  calculated  in  Figure  M-1 . Tints  of  the  radiance  as  function  of 
latitude  and  time  of  year  for  various  altitudes  is  shown  in  Figure  M-2.  The 

model  for  GPS  simulation  is  deduced  from  these  results.  For  maximum  radiance 
variation  a January  situation  is  considered.  The  radiance  is  assumed  to  be 
a function  of  latitude  but  not  longitude.  The  radiance  variation  for  south- 
ern hemisphere  is  assumed  to  be  that  for  northern  hemisphere  six  months 
later.  The  radiance  variation  for  earth  disk  is  approximated  by  piecewise 
linear  segments  using  the  values  corresponding  to  20  Km  altitude  shown  in 
Figure  M-2.  The  radiance  values  and  the  normalized  variations  are  summarized 
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Figure  M-2.  Normalized  Earth  Radiance  Variation 
with  Latitude  at  January 
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In  Table  M-1.  The  normalized  variation  is  used  to  be  compatible  with  the 
existing  earth  sensor  simulation  that  does  not  include  the  details  of 
sensor  optics  and  filter  response.  The  radiance  for  the  15°  latitude  band 
is  taken  as  that  corresponding  to  a 220°K  earth  which,  in  conjunction  with 
the  sensor  optics  and  filter  response,  provides  3.15  ww  when  the  "A"  detector 
FOV  is  fully  illuminated. 


TABLE  M-1.  EARTH  RADIANCE  VARIATION  WITH 
LATITUDE  AND  TIME  OF  YEAR 


LAT. 

MON. 

W/M^ 

STER 

NORMALIZED 

variation 

-75° 

JAN 

4. 13 

0.72 

1 

o^ 

o 

0 

JAN 

4.63 

0.80 

o 

1 

JAN 

4.88 

0.85 

-30° 

JAN 

5.5 

0.96 

+15° 

JAN/JUL 

5.75 

1.00 

30° 

JUL 

5.94 

1.03 

45° 

JUL 

6.19 

1.08 

60° 

JUL 

6 .44 

1.12 

75° 

JUL 

6.94 

1.21 
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Implementation  of  the  earth  radiance  variation  into  the  sensor  model  requires 
the  following  computations  which  are  based  on  the  geometry  given  in  Figure 
M-3. 


(a)  LOS-vector  for  detector  cell  element  centered  at  (X  , Y ) : 

3,  3. 


u 


2 2 
X + Y 
a a 


^s  = 


X sin  U 
a 


Y « 
s 


-Y  s i n U 
a 


Z = cos  U 
s 


Denote 


(b)  Look  point  on  earth  surface  along  L : 

a 


^^a  ^ ^IB  ^SB  ^-a  ^ 

(l/)  = (V)  + s (l/) 


whe  re 

(R  = look  point  on  earth  surface  along  L . 
“"3  a 

= vehicle  position  from  center  of  earth 
S = slant  range  from  ^ to  R 


3590.95  n.  mi. 
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I 

I 

i 

I 

I 

f 

1 


(c)  Look  point  latitude: 


Denote 


(I,-,  = 


R, 


1*^3 


S = sin  ^ 


look  point  latitude 


(d)  Scaling  factor  for  radiance  variation: 

V = F(E  ) 

where  F(6 ) is  the  piecewise  linear  function  defining  the 
radiance  variation  with  latitude;  F(ts)  is  plotted  in  Figure  M-1. 


(e)  The  value  of  AREA  is  modified  from  the  original  value  by  the 
scaling  factor,  i.e.: 

AREA-^AREA  * V 


Figure  M-4  presents  a comparison  of  the  GPS  simulation  radiance  model  and  a 
model  described  in  a Barnes  Engineering  document,  "Technical  Description 
Model  13-16901  and  Model  13-16903  Synchronous  Altitude  Static-Sensors", 

22  August,  1975.  Models  agree  reasonably  well. 
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BARNES'  MODEL 


Figure  M-4.  Comparison  o»  *’odels 
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SENSOR  ERROR  MODELING 

f 

The  sensor  error  modeling  considered  here  includes  detector  and  pre-amp 

noises,  electronics  parameter  shift  error,  error  due  to  temperature  ^ 

gradient  across  detector  and  sensor  alignment  error.  Horizon  uncertainty 

errors,  though  not  a sensor  originated  error,  are  also  included  in  the  con-  ^ 

sideration.  These  are  errors  that  exist  in  all  operational  situations. 

I 

Errors  due  to  sun  or  moon  in  sensor  detector  field  of  view  and  errors  | 

resulted  from  variation  of  earth  disk  size  due  to  vehicle  orbit  eccentricity 
exist  only  for  a specific  operational  condition  and  are  not  included  here. 

The  detector  and  preamplifier  noise  has  been  calculated  to  be  0.379  pV  (3c) 
for  sensor  in  track  mode  (A  detectors  only)  and  0.464  yV  (3c)  for  sensor 
in  acquisition  mode  (A+B  detectors)  with  a detector  scale  factor  of  3.15  pW/ 
deg.  and  responsivity  of  20  V /U . These  noises  are  equivalent  to  0.0060  deg. 

(3c)  and  0.0074  deg.  (3a)  respectively.  Use  the  output  scale  factors  of 
5 V/deg  and  1.25  V/deg  for  track  and  acquisition  modes  respectively.  The 
noise  at  output  then  becomes  0.03  V(3g)  and  0.0093  V(3o).  The  noise  error 
will  be  simulated  as  white  noise. 

i 

Electronics  parameter  shift  has  been  analyzed  by  Barnes,  Reference  M-2. 

The  values  at  initial  setting  and  at  end  of  life  have  been  calculated  to  be  f 

0.0034  deg.  (3o)  and  0.0219  deg.  (3o)  respectively.  These  correspond  to 

output  errors  of  0.017  V (3c)  and  0.1095  V (37)  for  track  mode  and  ' 

0.0043  V (3o)  and  0.0274  V (3a)  for  acquisition  mode.  The  error  due  to 

electronics  parameter  shift  will  be  simulated  as  bias  errors. 

The  temperature  gradient  considered  by  Barnes  amounts  to  a temperature 

o ^ 

difference  of  0.00048  F between  A and  S reference  junctions.  The  correspond- 
ing attitude  error  evaluated  for  temperature  of  /0°C  is  0.0103°  or  0.0515  V 
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for  track  mode.  For  acquisition  mode,  the  temfierature  change  of  0.0028  deg. 
between  detectors  is  calculated  to  cause  attitude  error  of  0.0029  deg.  which 
amounts  to  output  error  of  0.0035  V.  All  values  quoted  above  are  considered 
to  be  for  worst  case  and  will  be  treated  as  3°  values.  The  error  due  to 
tenperature  gradient  will  be  simulated  as  bias  errors. 

The  achievable  static  sensor  alignment  accuracy  has  been  estimated  by  Barnes 
to  be  in  the  order  of  0.015  deg.  (3a).  The  equivalent  output  error  is 
0.075  V (3c)  for  track  mode  and  0.019  V (3a)  for  acquisition  mode.  The 
effect  of  alignment  error  will  be  simulated  as  bias  errors.  The  horizon 
uncertainty  error  has  been  treated  as  an  exponentially  correlated  noise 
with  standard  deviation  0.88  Km,  correlation  time  10  days  and  correlation 
distance  2500  N . Mi  in  Reference  M-1,  i.e.. 


E[h(t,s)h(c+T,  S+d)]  = 


-d/2500 

e 


where 

a = 0.88  Km 
h 

T = time  in  days 
d = distance  in  n.  mi. 


— T /lO 

The  term  e can  be  approximated  by  1 if  T is  much  less  than  10  days. 

Factoring  in  the  GPS  orbit  altitude  of  10897  N.  Mi  the  horizon  uncertainty 
error  amounts  to  0.0020  deg.  of  earth  disk  angle  error.  In  the  evaluation 
of  Pj  and  detectors  at  opposite  sides  of  the  earth  disk  are  involved. 
The  separation  between  the  look  points  of  these  detectors  is: 


d = 


(180°  - lA°)X7r 
^o 


X 3461.6  N.  Mi. 


180 

= 100  29  N.  Ml. 
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The  separation  is  approximately  four  times  the  correlation  distance. 

Hence  the  effects  of  horizon  uncertainty  error  at  opposing  detectors 
are  essentially  uncorrelated.  The  equivalent  error  in  evaluation  of 
and  is  then  0.0020*^  x i'2  = 0.0028°  (lo)  or  0.0085°  On).  The  output 
error  becomes  0.0424  V for  track  mode  and  0.0106  V for  acquisition  mode. 

The  same  values  can  be  used  for  P2  and  R2 . (The  distance  between  look 
points  of  every  other  detector  cell  is  approximately  5383.  1 N .Mi . The 
correlation  of  horizon  uncertainty  at  points  separated  by  this  distance  is 
only  0.12.  It  is  neglected  for  simplicity  of  model.)  For  vehicle  points 
at  nadir  and  spins  at  a yaw  rate  of  W RPM,  the  sensor  detector  scans  the 
earth  surface  at  a rate  of  358.87  W N. Ml/sec.  The  250C  N.  Mi  correlation 
distance  would  then  be  equivalent  to  a correlation  time  of  6.97/W  sec. 

The  sensor  output  error  e(t)  due  to  the  horizon  uncertainty  will  be  modeled 
as  an  exponentially  correlated  noise  as  follows: 

E[e(t)  e(t+T)]  = e“^/’^c 

where 

Oq  = 0.0424  V for  track  mode 

= 0.0106  V for  acquisition  mode 

T = 6.97/W  sec  | 

c 

= 3. 485  sec  for  W = 2 RPM 

I 

The  above  model  will  be  used  to  account  for  the  effect  of  horizon  uncertainty 

even  for  cases  where  vehicle  is  not  quite  pointing  at  the  nadir.  (The  change  | 

of  correlation  time  due  to  the  pointing  offset  from  nadir  by  + 4°  is  expected 

to  be  small.  For  large  offset  from  nadir  the  effect  of  detector  not  seeing  i 

the  earth  disk  will  be  much  more  significant  than  tlie  effect  of  horizon 

uncertainty.) 

In  summary,  the  sensor  error  will  be  simulated  as  an  additive  output  error 

consists  of  a white  noise  component,  a bias  component  and  an  exponentially  ' 
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correlated  random  noise.  The  white  noise  is  used  to  simulate  the  effect 
of  detector  and  preamplifier  noises.  The  eqiiivalent  attitude  error  of  the 
white  noise  is  0.0060°  (3a).  The  exponentially  correlated  noise  is  used  to 
simulate  the  effect  of  horizon  xmcertainty  error.  The  equivalent  attitude 
error  is  0.0085°  (3a).  The  correlation  time  of  this  error  for  the  vehicle 
spin  rate  of  2 RPM  at  initiation  of  earth  acquisition  is  3.485  sec. 

The  bias  error  is  introduced  to  account  for  the  total  effects  of  electronics 
error,  tenperature  gradient  error  and  sensor  alignment  error.  The  equivalent 
attitude  error  is  0.0285°  (3o)  evaluated  as  the  RSS  value  of  eacli  individual 
e ffects . 

CONCERNS 

o The  detector  geometry  and  output  coiiT)utations  may  not  be  correct 
since  Barnes  has  tx^7o  different  designs.  The  one  used  here  (Tigure 
M-5)  Is  one  of  their  designs.  The  other  has  detector  pattern  rotated 
22^2°  about  the  line  of  siglit  from  the  one  shown  here. 

o The  output  update  frequency  for  GPS  is  not  certain.  The  one 
used  here  is  the  one  proposed  by  Barnes  for  DSCS  III. 

o Earth  image  distortion  and  iocus  are  not  modeled,  llxcludinj;  radiance, 
the  earth  image  is  assumed  to  be  a perfect  circle  on  the  focal  surface. 

o The  detector  output  approximation  using  1°  x 1°  squares  causes 
very  little  error  near  null  (normal  operation)  and  up  to  0.15° 
maximum  detector  error  at  extreme  attitudes  (during  acquisition). 

The  extent  of  the  errors  near  null  and  the  sensor  output  errors 
during  acquisition  has  not  been  analyzed.  The  approximation 
errors  off  null  are,  on  the  order  the  » • i ors  due  to  eartli  radiance 
variation  and  earth  image  qxidlity.  St  I ler  squares  may  be  used  to 
improve  accuracy  but  computer  cost  will  go  up. 
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NUMBERING  OF  X Y COORPI NATES  FOR  PORTION  OF  1°  \ T Si'UAR'^S. 
X-Y  COORDINATES  OF  ALL  160  f x l"  SQUARES  ARF.  03TAINE''  FROL. 
THESE  20  STORED  VALUES 


Figure  M-5.  Detector  eometvy— 2--Axis  i ih  Sensor  (cont  inin.  a') 
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Figure  M-5.  Detecloi 


metry— ?-/'x  is  Kartii  Sensor  (conc1u<Uvi) 


o Sun  and  moon  interference  are  not  included  in  this  model 


description. 

flow  diagram 

The  confcined  earth  sensor  model  flew  diagram  is  shown  in  Figure  M-6 . The 
left  arrow  within  boxes  indicates  replacement,  i.e.,  the  expression  on 
the  ri^t  replaces  the  variable  on  the  left.  Mathematical  subscripts  are 
used.  Symbols  used  in  the  flow  diagram  are  defined  in  the  sections  below 
Comments  appear  outside  the  boxes  to  describe  the  operations  within  the 
boxes . 

INITIALIZATION  DATA 

Table  M-2  lists  the  initialization  parameters,  nominal  values,  tolerance, 
and  source  of  data.  This  data  is  to  be  input  at  the  start  of  each  run. 

INPUT  DATA 

Table  M-3  lists  the  input  data  the  earth  sensor  model  expects  from  the 
calling  program  each  time  the  subroutine  is  executed  in  the  normal  mode. 

OUTPUT  DATA 

Table  M-A  lists  the  output  data  from  the  earth  sensor  model  subroutine. 
Other  subroutine  variables  may  be  outputs  if  the  user  desires. 

INTERNAL  VARIABLES 

^ Internal  variables  used  in  iiie  flow  diagram  art-  defined  in  Table  M-5. 

I 
I 
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SUBROUTINE 


Ini ti al ize 
time  counters 


Clear  detector 
outputs  and 
saved  detector 
area  coverages 


TYM  0 
TYMUP  0 


For  i 

= 1 thru  24 

Vi 

0 

DSAV^ 

0 

Ispin 
Mode) 
A Jp.225 


Input  values  for  X-j  thru  X^g 
Y]  thru  Y^o 
LOS^,  LOS^^, 

T^B  , At , SIG-j  thru  SIG-jg 
SFi  and  $Fo 
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p.224 


C0N2 


C0N3 


Detector  dynamic 
response  increments 


RA 

0, 

RIA 

0. 

R2A 

0 

RT 
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RIT 

0, 

R2T 

0 

PA 

0, 

PIA 

0, 

P2A 

0 

PT 

0. 

PIT 

0, 

P2T 

0 

Clear  outputs  of  2-axis  sensor 


STATE  1 


Set  initial  state 
of  2-axis  sensor 


AJ  p.225 


Figure  M-n.  Earth  Sensor  Flow  Chart 
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Figure  M-6.  ),arth  Sensor  Flow  i li.irt  (continued) 
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Figure  M-6 . ..arth  Sen.sor  Flow  i i.art  (continued) 


Figure  M-6 . liarth  Sensor  Flow  ( hart  (continued) 
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Figure  M-6,  Earth  Sensor  Flow  Htart  (continued) 
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Figure  M-ft.  Earth  Sensor  Flow 
Chart  (continued) 
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Figure  M-b,  i arth  Sensor  Flow  ( l.art  (continued) 
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Figure  M-6.  Earth  Sensor  Flow  Chart  (concluded) 
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Where 

SIGj  = 0.0020  (deg.,  1')  I = 1,  ....  8 

BSj  = Bias  error  with  worst  case  value  0.0285  (deg.,  3r)  1=1,  ...,  8. 

REXPj  = Exponent! al3.y  correlated  noise  with  standard  deviation 
0.00283°  and  correlation  time  3.A85  sec. 

The  function  REXP^  to  be  used  to  generate  the  exponentially  correlated  noise 
is  given  in  the  following: 

FUNCTION  REXPj 

REXPIO  = REXPl 

BETA  = 1 . /TMCNST 

DREXP  = -BErrA*REXPI0 

+SQRT  (2.*BETA)*RSIGI*RANDl 

REXPl  = REXPIO  + DREXP 

RETURN 

END 

Whe  re 

TMCNST  = Correlation  Time 
RSIGl  = Standard  Deviation 

RANDI  = Random  nuTiber  generator  with  unit  variance 
REXPl  to  be  assigned  an  initial  value. 
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TABLE  M-2.  INITIALIZATION  DATA  FOR  EARTH  SENSOR  MODEL 


PARAMtTtR 

SVBSCKl  PT 

— 

UNIT 

NOMINAL 

VALUE 

— 

SOURCE 

TOI,ERA.NCE  ■ SOURCE 

Tiine  constant  of 

two  axis  earth 

Seconds 

1.4 

1 ' t e I 

N/A  None 

sensor  detectors,  T, 

1 

Line  of  sight  of 

spin  earth  sensor 

field 

LOS°,  UXS 

1 



-0.984H 

Note  2 

N 'A  Ni*t 

■ 

speci f ied 

— 

-0.9848 

LOSy 

1 

0 

1 

0 

1 

mSj,. 

I 

0.  17> 

2 

-fM73A 

Trans  format i on  f ron 

hody  to  sensor 

frame:  (TSB) 

Row  1 Col.  1 

1 

— 

l.O 

AssuKieil 

L A '.eiie 

1 2 

’ 

0 

i nurdnaij 
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i 3 
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1 

2 1 
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2 2 

s 
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2 3 

n 

0 

1 

3 1 

0 

3 2 

n 

1 

3 3 

1 .0 
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Calling  t i 'X- 

‘ 

Increment,  t 

Set  onds 

'0.0032A 

in  m<’de  ' 

Sf*-  I 

text 

-0.125 

1.  axis) 

See 

text 

Center  pcsi tion 

of  Incremental 

dete<tor  are.i: 

*1 

I 

K/juiv. 

0.25307274)5  \.>t. 

N 1 ^.one 

2 

radi.in 

0.25  i07274r>  1 & 

1 

0..'7n52* 

Ill.l  ri>;.  3 

1 

I 


TABLE  M-2.  INITIALIZATION  DATA  FOR  EARTli  SENSOR  MODEL  (continued) 


ION  DATA  FOR  EAR'j 


I 


TABLE  M-2.  INITIAL! 


INSOR  MODEL  (concluded) 


Noise  Standard 

Deviation 

Sin  on:  RIT 

1 

Volts 

0.167 

Note  4 

N/A 

RlA 

2 

Volts 

0.0417 

PIT 

3 

Volts 

0.167 

PIA 

4 

Volts 

0.0417 

R2T 

5 

Volts 

0.167 

R2A 

6 

Volts 

0.0417 

P2T 

7 

Volts 

0.167 

P2A 

8 

Volts 

0.0417 

TPIP^ 

9 

Sec. 

0 

Not 

Sped  f- 
led 

Net 

Spec  1 ( 
led 

TPIP, 

10 

Sec. 

0 

Not 

Specif- 

ied 

Not 

Sped  f* 

it'd 

Two  axis  output 

scale  factor 

Track,  SF 

1 

Volt/sq . 
deg. 

1.25 

Note  3 

0.0625 

•Note  3 

Acq.  , SF 

2 

Volt/sq. 

0.3125 

0.025 

NOTES  TO  TABU  M-2 


Barnes  presentation  to  Honeywell  in  St.  Petersburg  on  3-l>*75. 
Rockwell  Spec  MC  432-0214,  Rev.  B.,  12-11-74, 

Computed  from  data  In  Rockwell  Spec  MC432-0214,  Rev.  B,  12-11-74. 
Track  scale  factor 


1 degree  attitude  ■ 4 sq.  deg.  coverage  * 5 l 0.25V. 

. 1.25  , 0.0625  volt/sq.  deg. 

4 sq . deg.  ^ * 

Acquisition  scale  factor 

4 degree  attitude  • 16  aq.  deg.  coverage 

‘--v’*  - - 0.3125  - 0.025  volt/sq.  deg. 

4 sq . deg.  ^ ^ 

Cotn>uted  from  data  In  Rockwell  Spec  MC  432-0214,  Rev.  B,  12-11-74 
T rack  mode 

tO. 1°,  3 sigma 

0.1  X 5 volt /deg  X 1/3  ■ 0. 16  7 volt  , 1 slpna 

Acquisition  mi  ile 

*•  . 1°,  3 sigma 

0.1  X 1.25  volt/deg  X 1/3  ■ 0.0417  volt,  1 sigma 
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TABLE  1^3.  INPUT  DATA  FOR  EARTH  SENSOR  MODEL 


definition 

1 

FLOW 

DIAGRAM  SYMBOL 

DIMENSION 

TYPE 

UNITS 

Mode  control  flag 

MODE 

Integer 

0 

= initialize 

1 

= spin  opera- 

tion 

2 

= 2 axis 

operation 

Satellite  position 

vector 

Transformation 

R ,R  ,R  y R 

X*  y 2 

3 

Real 

N 

Mi 

from  body  to 

TIB 

9 

Real 

inertial  frame 

L_ 
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TABLE  M-A.  i'  ' i’PUT  DATA  FOR  EAR  SEIiSOR  MODEL 


DEFINITION 

FLOW 

DIAGRAM  SYMBOL 

DIMENSION 

TYPE 

- 

UNITS 

Spin  earth  sensor 

output 

TP  IP 

2 

Real 

Seconds 

Pitch  attitude 

in  track  mode 

PT 

Real 

Volts 

Pitch  attitude 

in  acquisition 

mode 

PiV 

Real 

Volts 

Roll  attitude  in 

track  mode 

RT 

Real 

Volts 

Roll  attitude  in 

acquisition  mode 

RA 

Real 

Volts 
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TABLE  M-5.  INTERNAL  VARIABLES  FOR  EART  1 SENSOR  MODEL 


DEFINITION 

aow 

:iIAGR/i.M  SYMBOL 

DIMENSION 

TYPE 

UNITS 

Two  axis  rise/decay 

incretnenc  for  n.2'j  sec. 

CON  2 

Real 

Volts  per  sq . deg. 

Two  axis  rise/ decay 

CON  3 

Real 

Volfs  per  sq . deg. 

increment  for  0.125  sec 

Area  of  incremental 

A 

80 

Real 

Square  degree 

detector  covered  by 

earth 

Time  since  start  of  run 

TYM 

Doii) 

Seconds 

p rec 

Time  cl  next  2-axis 

lYMUP 

Doub 

Seconds 

update 

prec 

Satellite  distance  to 

R.M 

Real 

N.  Mi. 

earth  center 

Satellite  position 

vector  in  sensor 

frame 

RS  , RS  , RS  ; 

X y z 

R 

3 

Real 

.N.  Mi 

Radius  of  earth  image 

RE 

Real 

Equiv.  radian 

on  focal  surface 

Position  of  earth 

B 

Real 

I.quiv,  radian 

image  on  focal 

surface,  X rorponent 

Position  of  earth 

C 

Real 

Equiv.  radian 

Image  on  focal 

surface,  Y confionent 

State  of  seqiential 

processing  within 

STATE 

Intege  r 

sensor 

Area  of  <k  tector 

D 

12 

Real 

Square  degree 

co’.'ered  by  earth 

Area  of  detector 

coverage  from 

previous  cycle 

DSAV 

24 

Real 

Square  degree 

Dynamic  Detector 

V 

24 

Real 

Square  Degree 

output 
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TABLE  M-5.  INTERNAL  VARIABLES  FOP  EARTH  SENSOR  MODEL  (concluded) 


1 


DEFINITION 

FL(»W 

DIAHRAM  SYMBOL 

1 

DIMENSION 

TYPE 

UNITS 

Pitch  1 track 

PIT 

Real 

Volts 

attitude 

Pitch  2 track 

p:t 

Real 

Volts 

attitude 

Rcll  1 track 

RIT 

Real 

Volts 

attitude 

Roll  2 track 

R2T 

Real 

Volts 

attitude 

Pitch  1 acquisl- 

PlA 

Real 

Vol  ts 

tion  attitude 

Pitch  2 acqulsi- 

P2A 

Real 

Volts 

tlon  attitude 

Roll  1 acqulsi- 

RIA 

Real 

Volts 

tlon  attitude 

Roll  2 acqulsi- 

R2A 

Ri  al  V 

Volts 

tion  attitude 

Index 

1 

I n te  ge  r 

Index 

1 

Integer 

I.ine  of  sl^t  of  spin 

i.os' 

3 

Real 

sensor  field  in 

inertial  frar® 

Normalized  cross 

SANn 

Real 

product  maftnitude  of 

R and  LOS 

Komalized  dot  product 

GANG 

Real 

of  R and  1.0S 

Angle  between  -R  and 

A-VGI 

Real 

Radian 

tangent  to  earth 

Angle  between  -R  ;jnd 

AJir,2 

Real 

Radi  an 

r.os 

Angle  between  LOS 

/\NG3 

Real 

Radi  an 

and  earth  tangent 

Saved  angle  between 

ANGSAV 

2 

Real 

Radian 

earth  tangent  and  LOS 

Half  dimension  of 

X 1®  square 

RD 

Real 

Equlv  radian. 

(a  constant  preset 
to  0.0087266) 

Area  of  earth  cover- 

AREA 

Real 

Square  Degree 

Ing  1°  X 1°  square 

Funct  Ion 
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